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SUMMARY 

This  study  comprises  an  experimental  investigation  of  the  chemistry  of 
boron  combustion  products  with  particular  emphasis  on  the  phenomena  preceding 
and  accompanying  nucleation  in  expanding  flows.  Supersonic  (Mach  number  % 3) 
jets  of  the  burnt  gases  from  H2/O2/N2  flames  seeded  with  BC13  were  examined 
using  laser  light  scattering  and  a mass  spectrometer  equipped  with  a modulated 
molecular  beam  sampling  system.  The  results  demonstrate  that  the  processes  by 
which  the  primary  combustion  product  HB02  is  converted  to  condensed  phase  jet 
constituents  such  as  B203  (liquid)  are  not  interpretable  via  models  based  on 
conventional  homogeneous  nucleation  theories.  Mass  spectrometric  probing  of 
systems  held  close  to  conditions  favoring  particle  formation  reveal  that  the 
reaction  mechanisms  are  extremely  complex  and  involve  a wide  variety  of  inter- 
mediate species  of  widely  varying  H/B/O  proportions.  The  dehydration  corre- 
sponding to  the  overall  process  2HB02  -*•  B203  (liquid)  + H20  is  found  not  to  occur 
at  any  discreet  step  in  the  reaction  scheme,  but  rather  appears  to  occur 
gradually  and  is  still  proceeding  even  when  the  bulk  of  the  boron  present  is 
in  the  form  of  species  with  molecular  weights  between  150  and  500  amu.  The 
mechanistic  path  to  molecules  with  molecular  weights  in  excess  of  135  seems  to 
proceed  through  addition  of  the  'monomers'  HB02,  B203  and,  less  often,  H3B03 
to  a series  of  substrates  of  which  H2B305  is  the  lowest  member. 

Rate  coefficients  for  the  reactions  leading  to  the  dimer  and  trimer  of 
HB02 , vJ  t. , 

2HB02  -►  (HB02)2  and  HB02  + (HB02)z  -*•  (HB02)3 

have  been  determined  at  a temperature  of  650  K to  be  8 x 10-11  and  3 x IO-10 
ml  molecule-1  sec-1,  respectively.  The  HB02  polymerization  process  appears 
to  terminate  with  the  trimer,  however,  and  neither  it  nor  the  dimer  appear  to 
play  an  important  role  in  the  overall  condensation  process. 

During  the  course  of  this  work  some  additional  experiments  on  dissociative 
electron  attachment  to  HB02  and  HC1  were  also  performed  to  augment  results  ob- 
tained during  prior  ONR  work.  The  data  have  been  interpreted  to  yield  a rate 
coefficient  of  3.2  x 10-l°  exp (-11 ,000/T)  ml  molecule-1  sec-1  for  the  reaction 
e-  + HBOj  -*•  B02-  + H and  a value  of  6 x 10-J1  exp (-10 , 100/T)  ml  molecule-1 
sec-1  for  e-  + HC1  -+  Cl-  + H. 
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I.  INTRODUCTION 


The  primary  goal  of  this  research  is  the  definition  of  the  reaction 
mechanisms  which  link  the  large  number  of  known  boron  combustion  products, 
including  condensed  phases,  to  one  another.  During  the  early  stages  of  the 
work,  however,  while  the  required  instrumentation  was  being  assembled  and 
modified,  additional  experiments  were  performed  to  augment  a body  of  data  on 
electron  attachment  to  HB02  and  HC1  which  had  been  collected  during  earlier 
ONR  and  SAMSO  work  (Contracts  N00014-74-C-0326  and  F04701-72-C-0079,  respec- 
tively). The  results  of  these  efforts  are  only  summarized  in  the  body  of  this 
report;  details  are  available  in  the  resultant  publications  included  as 
Appendices  A1  and  B.a  The  reprint  of  an  invited  review  paper  on  charged  spe- 
cies diagnostics  in  combustion  systems  prepared  under  the  auspices  of  this 
contract  for  the  14th  Aerospace  Sciences  Meeting  of  the  AIAA  is  also  included 
as  Appendix  C.3 


1.  Miller,  W.J.  and  Gould,  R.K.,  "Electron  Attachment  Kinetics  in  Flames: 
Dissociative  Attachment  to  HBO2,"  Chem.  Phys.  Lett.  38,  237-241  (1976). 

2.  Miller,  W.J.  and  Gould,  R.K.,  "Electron  Attachment  Kinetics  in  Flames: 
II.  Dissociative  Attachment  to  HC1,”  AeroChera  TP-337,  February  1976, 
submitted  to  J.  Chem.  Phys. 

3.  Miller,  W.J.,  "Charged  Species  Diagnostics  for  Combustion  Systems," 
AeroChem  TP-331,  AIAA  Paper  76-135,  presented  at  14th  AIAA  Aerospace 
Sciences  Meeting,  Washington , DC,  January  1376. 
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A.  DISSOCIATIVE  ATTACHMENT  TO  HBO? 

Addition  of  boron  compounds  to  high  temperature  combustion  flames 
containing  electrons  gives  rise  to  the  formation  of  the  negative  ion  B02“. 

The  thermodynamic  properties  of  this  ion  were  established4* ' 5 in  our  prior  ONR- 
supported  work  and  confirmed  in  subsequent  studies  by  others.6  Its  electron 
affinity  of  410  kJ  mole-1  (4.25  eV)  renders  it  an  attractive  potential  candidate 
for  the  control  of  electron  densities  in  wakes  and  exhaust  plumes  as  well  as  a 
considerable  possible  nuisance  in  MHD  ducts.  To  quantitatively  assess  the 
effects  of  its  formation  in  dynamic  systems,  however,  the  rate  of  its  formation 
must  be  known  in  addition  to  its  equilibrium  properties.  The  rate  coefficient 
for  the  dissociative  attachment 

e + HB0a  -*•  B02“  + H (!) 

was  therefore  measured  in  a number  of  standard,  well  characterized,  laboratory 
test  flames  at  total  pressures  of  0.13  and  1.0  atm.  The  data  are  consistent 
with  the  Arrhenius  expression  k = 3.2  x 10~10  exp (— 11 , 000/T)ml  molecule  1 sec  1 
between  1730  and  2250  K.  The  failure  of  the  observed  B02“  formation  rates  to 
correlate  in  any  systematic  way  with  hypothetical  three-body  rate  coefficients 
demonstrates  that  three-body  attachment  processes  are  negligible  in  these 
flames  and  that  dissociative  attachment  dominates. 

A more  detailed  account  of  these  measurements  may  be  obtained  in 
Ref.  1,  a reprint  of  which  is  included  as  Appendix  A. 


4.  Jensen,  D.E. , "Electron  Attachment  and  Compound  Formation  in  Flames. 

I.  Electron  Affinity  of  B02  and  Heats  of  Formation  of  Alkali  Metal 
Metaborates,"  Trans.  Faraday  Soc.  6_5,  2123-2132  (1969). 

5.  Jensen,  D.E.,  "Electron  Attachment  and  Compound  Formation  in  Flames. 

II.  Mass  Spectrometry  of  Boron-Containing  Flames,"  J.  Chem.  Phys.  52^, 
3305-3306  (1970). 

6.  Srivastava,  R.D.,  Uy,  O.M.,  and  Farber,  M. , "Effusion-Mass  Spectrometric 
Study  of  the  Thermodynamic  Properties  of  B0“  and  B02“,"  Trans.  Faraday 
Soc.  67,  2941-2944  (1971). 
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B.  DISSOCIATIVE  ATTACHMENT  TO  HC1 


The  highly  electrophilic  nature  of  chlorine  coupled  with  its 
ubiquity  render  its  negative  ion  chemistry  the  most  important  plasma  processes 
in  a wide  variety  of  high  temperature  environments.  Despite  the  obvious  signi- 
ficance of  Cl-  based  on  its  thermochemical  properties  and  high  concentration 
in  practical  systems,  there  was,  until  this  work,  no  direct  measure  available 
of  its  formation  rate  at  high  temperatures.  Results  of  some  of  our  prior  work 
were  therefore  re-analyzed  and  supplementary  experiments  performed  to  deter- 
mine the  mechanism  for  Cl-  formation  and  the  rate  coefficient  for  the  rate 
limiting  step. 

The  experiments  involve  mass  spectrometric  measurements  of  negative 
ion  concentrations  as  functions  of  flow  time  in  the  burnt  gases  of  well 
characterized  laboratory  test  flames  seeded  with  small  quantities  of  alkali 
metals  and  chlorine  compounds.  The  results  are  consistent  with  the  two-body 
dissociative  attachment  mechanism 

e + HC1  ->  Cl"  + H (2) 

with  a rate  coefficient  well  represented  by  the  Arrhenius  expression 

k = 6 x 10~11  exp(— 10,000/T)ml  molecule-1  sec-1  (3) 

over  the  temperature  range  1720  K to  2475  K.  The  flames  are  operated  at 
pressures  of  100  and  760  Torr.  Electron  concentrations  were  measured  using 
a microwave  resonance  cavity  method  and  positive  ion  profiles  (total  charged 
species)  were  determined  with  electrostatic  probes.  H-atom  concentration 
profiles  required  to  take  account  of  the  reverse  detachment  process  are  well 
known  for  the  760  Torr  flames;  in  the  low  pressure  flames  they  were  obtained 
using  a mass  spectrometric  negative  ion  tracer  technique  involving  the  known 
ion-molecule  chemistry  of  molybdenum  additives  and  confirmed  via  optical 
absorption  measurements  of  [OH].  Contributions  to  Cl-  formation  from  three- 
body  attachment  processes  are  not  observed.  The  measured  values  of  k are  found 
to  be  somewhat  lower  than,  but  in  general  agreement  with,  those  inferred  from 
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electron  beam7’8  and  swarm  studies9  but  are  much  lower  than  the  values  pre- 
viously derived  from  measurements10*11  of  the  reverse  reaction  rate. 


The  detailed  description  of  this  work  is  presented  in  Appendix  B. 
The  manuscript  is  currently  being  revised  somewhat  and  is  expected  to  be 
accepted  soon  for  publication  in  the  Journal  of  Chemical  Physics. 
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7.  Buchel ' nikova , I.S.,  "Cross  Sections  for  the  Capture  of  Slow  Electrons 
by  02  and  H20  Molecules  and  Molecules  of  Halogen  Compounds,"  Sov.  Phys. 
JETP  35,  783-791  (1959). 


8.  Azria,  R. , Roussier,  L. , Paineau,  R. , andTronc,M.,  "Attachement  Elec- 
tronique  Dissociatif  sur  HCl  et  DCi,"  Rev.  Phys.  Appl.  _9,  469-473  (1974). 

9.  Christophorou,  L.G.,  Compton,  R.N.,  and  Dickson,  N.W. , "Dissociative 
Electron  Attachment  to  Hydrogen  Halides  and  their  Deuterated  Analogs," 

J.  Chem.  Phys.  48,  1949-1955  (1968). 

10.  Howard,  C. J. , Fehsenfeld,  F.C.,  and  McFarland,  M. , "Negative  Ion-Molecule 
Reactions  with  Atomic  Hydrogen  in  the  Gas  Phase  at  296 °K,"  J.  Chem.  Phys. 
60,  5086-5089  (1974). 

11.  Burdett,  N.A.  and  Hayhurst,  A.N.,  "Kinetics  of  Formation  of  Chloride  Ions 
in  Atmospheric  Pressure  Flames  by  Way  of  HCl  + e~  t H + Cl“,"  Nature 
Phvr. . Sci.  245,  77-78  (1973). 
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III.  BORON  COMBUSTION  PRODUCT  CHEMISTRY 
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k.  PRACTICAL  CONSIDERATIONS 

A good  deal  of  work  has  been  performed  in  the  last  decade  on  the 
thermochemistry  of  boron  oxides  and  oxyhydroxides  and  on  the  overall  kinetics 
of  boron  metal  particle  combustion  (e.g.  see  Refs.  12-14).  The  stimulus  for 
these  efforts  is  the  very  attractive  theoretical  performance1*  of  boron  as  a 
high  energy  propellant  constituent.  One  critical  aspect  of  our  current  knowl- 
edge of  boron  combustion  chemistry,  namely  the  kinetics  and  mechanisms  of 
oxide  and  oxyhydroxide  interconversion  processes,  remains  poorly  understood, 
however.  It  has  been  demonstrated11*  that  the  presence  of  water  vapor  above 
a burning  boron  particle  greatly  enhances  product  volatility  and  burning  rate 
by  converting  liquid  B203  on  the  surface  to  gaseous  HB02.  The  interconversion 
of  these  two  combustion  products  is  also  involved  in  important  heat  release 
processes  which  account  for  a significant  portion  of  the  potentially  attain- 
able propulsion  performance.  As  noted  above,  the  objective  of  the  bulk  of 
the  research  described  in  this  report  is  the  elucidation  of  the  chemistry  of 
these  boron  combustion  products  through  an  investigation  of  the  mechanisms  by 
which  HB02  reacts  to  produce  condensed  phase  products  in  propulsion  system 
exhaust  streams;  the  oxidation  of  the  metal  itself  is  not  addressed  in  this 
work.  However,  it  should  be  noted  that  many  of  the  processes  involved  in 
these  product  interconversion  mechanisms,  or  their  reverse  reactions,  must 
also  comprise  important  steps  in  metal  combustion  at  higher  temperatures. 


12.  Macek,  A.  and  Semple,  J.M. , "Combustion  of  Boron  Particles  at  Elevated  3 

Pressures,"  Thirteenth  Symposium  (International)  on  Combustion  (The  Com-  i 

bustion  Institute,  Pittsburgh  1971),  pp.  859-868. 

13.  Macek,  A.  and  Semple,  J.M. , "Combustion  of  Boron  Particles  at  Atmospheric 
Pressure,"  Combust.  Sci.  Techn.  1_,  181-191  (1969). 

14.  Roberts,  R. , "Boron  Combustion:  A Review,"  Proceedings  of  the  9th  Inter- 

national Symposium  on  Space  Technology  and  Science,  Tokyo  1971,  pp.  123-135. 

15.  All  thermodynamic  data  taken  from  JANAF  Thermochemical  Tables,  Dow  Chemi- 
cal Co.,  Midland,  Mich.  ' 
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The  energetic  implications  of  achieving  the  formation  of  condensed 
phase  products,  namely  B203  in  the  liquid  state,  may  be  demonstrated  by  con- 
sidering the  thermodynamics  of  the  alternative  overall  metallic  boron  oxida- 
tion schemes 


2B  + 202  + H2 


■B203 (liq)  + H20 


AH20oo  = -275.4  kcal  (4) 

AH200o  = -351.2  kcal  (5) 


At  a representative  propulsion  system  exhaust  temperature  of  2000  K (that 
temperature  for  which  these  heat  changes  are  given)  the  equilibrium  constants 
(Kp)  for  these  two  processes  are  nearly  identical  and  both  products  are  there- 
fore of  nearly  equal  importance.*  If  condensation  fails  to  occur  and  B203(liq) 
is  not  formed,  roughly  half  of  the  76  kcal  difference  in  these  two  pathways 
is  not  realized  and  the  system  suffers  approximately  a 10%  decrease  in  over- 
all performance.  At  lower  exhaust  temperatures  the  equilibrium  concentrations 
of  B203 (liq)  are  even  higher  with  respect  to  HB02  and  if  condensation  still 
does  not  occur,  the  loss  of  performance  is  correspondingly  greater. 

To  further  quantify  these  effects  and  determine  their  significance 
in  practical  systems,  we  have  computed  the  specific  impulse  (Igp),  with  and 
without  B203  condensation,  for  several  air  augmented  solid16  and  air  breath- 
ing liquid17  propellant  systems  under  consideration18  in  the  late  1960's. 


* This  is  true  only  to  a first  approximation,  of  course;  the  exact  ratio  of 
these  and  other  B-containing  products  at  equilibrium  depends  on  the  detailed 
stoichiometry,  B content,  pressure,  etc.  of  the  entire  system. 

16.  Deklau,  B. , et  al,  "Design  Considerations  for  Air-Augmented  Rocket  Missiles 
and  Related  Afterburning  Testing"  (U) , 3rd  ICRPG/AIAA  Solid  Propulsion 
Conference  (U) , CPIA  Publication  No.  167,  Vol.  I,  April  1968,  pp,  481-500. 
(CONFIDENTIAL) 

17.  Rosenberg,  S.D.,  Yates,  R.E.,  and  Adrian,  R.C.,  "Secondary  Combustion  of 
Pentaborane  - Hydrazine  Exhaust  in  Air"  (u)  , 10th  Liquid  Fropulsion 
Symposium  (U) , CPIA  Publication  No.  176,  Vol.  I,  October  1968,  pp.  565-577. 
(CONFIDENTIAL) 

18.  Bubb,  J.E.,  "Future  Requirements — Airbreathing  (U),"  Propulsion  Require- 
ments for  the  70' s,  AIAA  6th  Propulsion  Joint  Specialist  Conference,  CPIA 
Publication  No.  197,  July  1970,  pp.  13-58.  (CONFIDENTIAL) 
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The  conditions  chosen  for  the  performance  calculations  were  flight  Mach  num- 
bers 3.6  and  3.0  at  20,000  ft.  altitude  and  a nozzle  exhaust  pressure  of 
7 psig  (roughly  ambient).  The  calculations  were  run  with  the  Air  Force 
Rocket  Propulsion  Laboratory  (AFRPL)  ISP  Code  which  includes  an  air- 
augumentation  routine.  The  results  are  summarized  in  Tables  1 and  2. 

TABLE  1 

AIR  AUGMENTED  SOLID  PROPELLANT 


Propellant : 

Air/Propellant  (wt) : 

Chamber  Pressure  (psia) : 

Flight  Mach  No. : 

Icri^sec^  with  B203  Condensation: 
(sec ) 

lSp  without  Condensation: 


Bi  , o C o . i»  O 3OI  o • 4 . OO3H0  • 735N0.OA2O0.  17  0 


Flight  Mach  No.: 

Isp^sec^  with  Ba03  Condensation: 
ISp(8ec)  without  Condensation: 


TABLE  2 

LIQUID  B0RANE/AIR  PROPELLANT 


Propellant: 

Air /Propellant  (wt): 
Chamber  Pressure  (psia): 
Flight  Mach  No.: 


(sec) 


with  B203  Condensation: 


ISp(sec)  without  Condensation: 


Bt . oH3 . s iN0 . 5 9 


Flight  Mach  No.: 

ISp(sec)  with  B203  Condensation: 
ISp^sec)  without  Condensation: 
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A great  many  of  these  propellants,  when  tested,  failed  to  achieve 
the  95-plus%  efficiencies  necessary  for  acceptable  performance.  It  was 
found16’ 17  that  large  air-propellant  ratios  were  required  to  get  good  per- 
formance and  the  limitation  was  tacitly  assumed  to  be  tied  to  incomplete 
metal  combustion.  Analyses  of  the  exhaust  products,  in  the  few  cases  where 
they  were  done,17  indicate,  however,  that  very  little  B203  (condensed)  is 
present  in  the  products.  This  strongly  suggests  that  at  least  a portion  of 
the  poor  performance  is  attributable  to  the  product  interconversion  chemistry 
which  is  the  subject  of  the  present  study. 

The  data  of  Tables  1 and  2 indicate  decreases  in  system  efficiencies 
of  about  7 and  5%  for  the  solid  and  liquid  propellants,  respectively.  Typical 
experimental  results  produce  deficiencies  in  the  range  10%-15%  or  more  at 
high  boron  loading  levels.  It  thus  appears  that  if  the  origin  of  the  condensa- 
tion losses  could  be  identified  and  overcome,  a substantial  fraction  of  these 
inefficiencies  could  be  eliminated. 

B.  APPARATUS 

In  these  experiments,  the  burned  gases  from  a number  of  H2/02/N2 
flames  containing  BC13  additive  are  expanded  through  a nozzle  and  the  charac- 
teristics of  the  supersonic  jet  emerging  from  the  nozzle  exit  are  examined  as 
a function  of  added  BC13.  The  use  of  a number  of  flames  allows  for  variations 
in  reservoir/jet  temperatures  and  bulk  gas  chemical  composition.  Our  effort 
comprised  two  sets  of  diagnostic  procedures;  (i)  laser  light  scattering  studies 
designed  to  define  the  limiting  conditions  for  condensation  and  (ii)  molecular 
beam  mass  spectrometry  in  boron-containing  combustion  product  jets  under  con- 
ditions up  to  and  including  those  defined  in  (i)  above. 

1.  The  Low  Pressure  Flame  Reservoir  and  Supersonic  Jet 

The  premixed  flames  which  provide  the  reservoir  section  for  the 
nozzle  expansion  are  members  of  an  extensive  set  of  subatmospheric  pressure, 
laminar  H2/02/N2  flames  used2  in  previous  kinetic  studies.  They  are  supported 
in  a Pyrex  vacuum  housing  on  an  annular  burner  with  an  inner  test  flame  matrix 
of  2.6  cm  diam  and  an  overall  diameter  of  A. 5 cm.  The  unburned  gas  flow  to 
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the  burner  is  20  to  30  ml (STP) sec-1  cm-2  (burner  area).  These  flows  provide 
stable  flames  in  the  50-200  Torr  pressure  range.  The  unburned  gas  composi- 
tions of  the  test  flame  series,  the  individual  temperatures  as  measured  by 
Na  D-line  reversal  and  the  computed  burnt  gas  equilibrium  compositions  of 
each  at  the  corresponding  measured  temperature  are  given  in  Table  3. 

The  BC13  additive  is  introduced  into  the  center  flame  N2  supply 
line  from  an  inverted  buret  flow  metering  system.  Ideally,  gas  flows  are 
tailored  so  that  only  the  central  test  flame  combustion  products  enter  the 
nozzle;  the  excess  gas  from  the  annular  shield  flame  is  pumped  off  through 
a bypass.  This  condition  is  not  always  attainable,  however,  and  some  dilution 
of  the  test  gas  by  a portion  of  the  shroud  must  be  accepted  when  using  lower 
velocity  flames. 

The  nozzle  through  which  the  test  flame  combustion  products 
are  expanded  is  contoured  to  give  a parallel,  shock-free  jet  with  a Mach  num- 
ber (M)  of  % 3.  Its  exit-to-throat  area  ratio  is  4.4,  the  throat  area  and 
throat-to-exit  distance  being  0.32  cm2  and  1.6  cm,  respectively.  Typical 
ratios  of  chamber  pressure  to  jet  static  pressure  under  conditions  of  perfect 
expansion  are  about  27  and  the  jet  gas  velocity  is  % 1.6  x 103  cm  sec-1. 


2.  Laser  Light  Scattering 

The  presence  of  particulate  material  in  the  jet  is  monitored 
optically19  by  observing  the  light  scattered  by  90°  from  an  incident  laser 
beam.  The  optical  apparatus  configuration  is  shown  in  Fig.  1.  The  jet  of 
combustion  product  gas  issuing  from  the  exit  of  the  supersonic  (Mach  3)  noz- 
zle is  intersected  by  a chopped  laser  beam  and  the  scattered  light  detected 
by  an  HTV  R212  photomultiplier  tube  (PMT)  focussed  on  the  point  of  laser 
beam/jet  intersection.  The  laser  light  at  732.8  nm  is  provided  by  a Spectra- 
Physics  Model  133  He-Ne  laser  and  is  modulated  at  1 kHz  with  a slotted  disc 
rotating  chopper.  The  PMT  output  signal  is  processed  and  displayed  by  a PAR 
Model  128  lock-in  amplifier. 


4 


] 


19.  Stein,  G.D.  and  Wegener,  P.P.,  "Experiments  on  the  Number  of  Particles 
Formed  by  Homogeneous  Nucleation  in  the  Vapor  Phase,"  J.  Chera.  Phys.  46, 
3685-3686  (1967). 
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Fig.  1 APPARATUS  FOR  PARTICLE  DETECTION  VIA  LASER  LIGHT  SCATTERING 


Experimentally,  the  optical  system  was  arranged  geometrically 
and  tested  in  the  absence  of  a flame  by  entraining  tobacco  smoke  in  the 
burner  chamber  supply  gas  stream.  The  capability  of  the  system  to  detect 
particles  in  actual  combustion  product  mixtures  was  then  established  by  seed- 
ing of  Hj/Oj/Na  or  Ar  flames  with  Fe(C0)5.  The  burner  housing  in  these  tests 
was  maintained  at  about  50  Torr  and  the  static  jet  pressure  is  1.5  ± 0.5  Torr. 

The  onset  of  condensation  was  readily  observed  as  a distinct 
increase  in  the  scattered  light  signal.  The  condensed  phase  material  detect- 
ed in  this  instance  however  is  probably  formed  in  the  flame  itself  rather 
than  in  the  nozzle  during  expansion;  equilibrium  calculations  of  the  burnt 
gas  composition  indicate  that  the  iron  should  be  present  in  these  flames  ■ 

largely  as  FeO(liquid)  or  Fe30<,  (.solid) . Such  is  not  the  case  for  boron  com- 
bustion products — thermodynamic  considerations  preclude  the  formation  of  \ 
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condensed  phase  boron  oxides  or  acids  anywhere  but  in  the  cooled  gases  down- 
stream of  the  nozzle  throat. 

3.  Molecular  Beam  Mass  Spectrometer 

The  molecular  beam/mass  spectrometer  apparatus  employed  in 
these  studies  was  designed  following  what  have  become  conventional  procedures 
with  the  exception  of  the  initial  sampling  stage.  Since  the  gas  to  be  sam- 
pled is  present  as  a supersonic  stream,  the  first  inlet  orifice  must  be 
designed  to  attach  or  swallow  the  shock  created  as  the  jet  impinges  upon  it 
and  thus  allow  for  undisturbed  sampling  of  the  jet.  For  the  inviscid  contin- 
uum supersonic  flow  over  a cone  at  zero  angle  of  attack  it  is  well  known20 
that  the  shock  wave  enveloping  the  cone  will  be  conical  and  attached  to  the 
cone  tip  only  if  the  cone  semi-angle  does  not  exceed  a critical  value,  0cri£» 
which  is  a function  only  of  the  free  stream  Mach  number  and  the  specific  heat 
ratio,  y,  of  the  incident  gas.  For  blunter  cones  the  shock  wave  becomes 
detached  (see  Fig.  2a).  Numerical  values  of  0cr£t  have  been  computed  for 
Y ~ l.A  with  the  results  for  2 Mach  numbers  given  in  the  first  row  of  Table  4. 
Referring  to  these  entries,  it  is  clear  that  if  we  are  interested  in  sampling 
in  the  Mach  number  range  2-3,  the  attached  conical  shock  condition  would  only 
be  met  if  the  overall  cone  angle  were  less  than  about  40°.  However,  a sampling 

TAF.LE  4 

VALUES3  OF  9crlt  FOR  yc-  = 1.4 


Geometry 


Moo  = 2 Me  = 3 


Cone  (3)  40.6°  49.1° 

Wedge  (2)  23.0°  34.0° 

g 

Estimated  from  Figs.  19  and 
22  of  Ref.  21. 


20.  Shapiro,  A.H.,  The  Dynamics  and  Thermodynamics  of  Compressible  Fluid  Flow 
(Ronald  Press,  New  York,  1953),  Vol.  1,  Section  16.5,  pp.  544-551. 

21.  Van  Driest,  E.R. , "Basic  Relations  in  Gas  Dynamics,"  Handbook  of  Engineer- 
ing Fundamentals,  O.W.  Eshbach,  ed.  (John  Wiley,  New  York,  1953)  2nd  Ed. 
Section  7,  Part  7,  pp.  7-08  - 7-39. 
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Fig.  2 SAMPLING  CONE  DESIGN 

a.  Cone/Shock  Interactions 

b.  Experimental  Configuration 
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cone  must  contain  an  orifice  and,  if  we  focus  our  attention  in  the  immediate 
vicinity  of  the  orifice,  the  flow  more  closely  resembles  (becomes  asymptoti- 
cally equivalent  to)  the  flow  over  a two-dimensional  wedge.  Interestingly 
enough,  supersonic  flow  over  a wedge  is  qualitatively  similar  to  that  over  a 
cone  in  that  a critical  wedge  angle  exists  above  which  the  flow  will  become 
detached.  If  one  examines  the  computed  values  of  9crit  f°r  t^ie  supersonic 
flow  over  a wedge  (see  row  2 of  Table  4)  it  is  found  that  they  are  signifi- 
cantly smaller  at  the  same  Mach  number.  Thus,  again  referring  to  Table  4,  it 
is  anticipated  that  the  required  cone  angle  in  the  immediate  vicinity  of  the 
orifice  would  have  to  be  kept  below  23°  in  the  Mach  2-3  range,  but  'far'  down- 
stream from  the  orifice  a cone  angle  of  40°  would  suffice.  We  are  thus  led 
to  a configuration  sketched  in  Fig.  2b  which,  in  the  absence  of  viscous  flow 
effects,  satisfies  all  of  the  attached  shock  conditions. 

Prior  work  on  ion  sampling  using  cones  of  the  above  design  has 
demonstrated22  that,  in  practice,  this  ideal  behavior  is  not  achieved.  How- 
ever, the  shock  present  at,  or  in  the  vicinity  of,  the  sampling  orifice  has 
been  found  to  be  weak  and  to  have  little  or  no  effect  on  the  sample;  this  is 
apparently  due  to  the  low  gas  densities  produced  as  the  sample  is  re- 
accelerated from  the  region  directly  in  front  of  the  orifice. 

The  skimmer  which  excludes  all  but  the  central  portion  of  the 
sampled  gas  stream  is  located,  as  shown  in  Fig.  3,  between  100  and  200  cone 
orifice  diameters  downstream  of  the  sampling  inlet.  The  outside  angle  of  the 
skimmer  cone  is  50°;  the  inside  angle  is  45°.  The  walls  terminate  at  the 
cone  tip  in  knife  edges  around  a 0.25  mm  diameter  hole.  Approximately  5 cm 
downstream  of  the  skimmer  opening,  the  beam  is  chopped  with  a normally  closed 
tuning  fork  chopper  (and  collimated  into  an  Extranuclear  Model  041-1  co-axial 
electron  impact  ion  source) . 


; 22.  Burke,  R.R.  and  Miller,  W.J.,  "Study  of  Mass  Spectrometric  Ion  Sampling 

1 Processes,"  Final  Report,  AeroChem  TP-247,  AFCRL-70-0550,  DDC  AD  725  149, 

j September  1970. 
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Fig.  3 MOLECULAR  BEAM  MASS  SPECTROMETER 


Following  ionization,  the  sample  ions  are  focussed  and  injected 
into  the  mass  filter.  The  signal  produced  at  the  ion  multiplier  tube  then 
feeds  a lock-in  amplifier  and  the  signal  is  displayed  on  an  oscilloscope  or  a 
strip  chart  recorder.  The  molecules  originating  in  the  sampling  system  are 
distinguishable  as  a 140  Hz  a.c.  signal  superimposed  (at  most  mass  peaks)  on 
a d.c.  component  resulting  from  the  ionization  of  background  gas  in  the  ionizer. 
The  a.c.  sensitivity  at  any  specific  mass  peak  is  dependent  on  the  magnitude 
of  this  background  and  the  identity  of  the  gas  but  it  typically  ranges  from 
about  10  parts  per  million  (ppm)  when  no  background  is  present  to  about  100  ppm 
where  a large  d.c.  signal  is  present. 
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1.  Laser  Light  Scattering  Condensation  Studies 


The  scattered  light  signal  from  the  laser  illuminated  jet  was 
determined  as  a function  of  BC13  loading  in  each  of  five  of  the  test  flames  of 
Table  3;  three  02-rich  flames  and  two  H2-rich  analogues  were  utilized  to  pro- 
vide a range  of  temperatures  and  variable  gas  composition.  The  two  coolest 
flames,  P6  (H2-rich,  n.  1550  K)  and  0P6  (02-rich,  'v  1500  K)  gave  clear  indica- 
tions of  particle  formation  at  about  1%  added  BC13  in  the  test  flame  gas  supply. 
In  the  three  hotter  flames,  P3  (H2-rich,  1730  K)  0P3  (02-rich,  1700  K)  and  0U4 
(02-rich,  n,  1650  K) , up  to  9%  BC13  was  added  without  observing  any  significant 
scattered  light. 


i 

1 

1 

i 

1 

1 

t 


\ 

j 

j 


Figure  4 shows  some  of  the  results  of  the  light-scattering  j 

experiments  conducted  using  the  nozzle  described  in  Section  III.B.  A burner  j 

housing  pressure  of  100  (±  10)  Torr  and  a static  pressure  of  2.0  (±  0.3)  Torr  \ 

were  used  to  insure  overexpansion  and  thus  minimize  the  possibility  of  inter-  j 

ference  from  shock  formation  in  the  nozzle  or  along  the  jet  axis  between  the  j 

exit  plane  and  the  laser /jet  intersection  point.  j 

j 

During  the  course  of  our  early  experiments,  us-'ng  a larger,  I 

0.64  cm  diam  nozzle,  under  conditions  of  high  BC13  addition  rates,  it  was  noted  j 

that  deposits  of  white  powdery  condensate  built  up  on  the  converging  section  of  | 

the  nozzle  inside  the  flame  reservoir  section  and  that  particles  were  present 
at  the  cool  outer  boundaries  of  the  shroud  flame. 

It  was  thus  not  certain  whether  all  of  the  observed  condensed 
phase  was  formed  in  the  jet  or  whether  some  of  it  was  formed  within  the  reser- 
voir housing.  To  remove  this  ambiguity,  the  0.64  cm  diam  nozzle  was  replaced 
with  one  with  a throat  diam  of  0.32  cm.  The  smaller  nozzle  (i)  allows  opera-  ; 

tion  at  higher  burner  pressures  (which  is  advantageous  from  the  point  of  view  \ 

of  promoting  condensation)  and  (ii)  allows  a large  amount  of  the  burned  gases  3 

to  be  pumped  away  through  a bypass  valve  while  only  gases  from  the  core  of  the 
flame  pass  through  the  nozzle.  These  two  nozzles,  although  they  necessitate 
the  use  of  quite  different  reservoir  and  jet  pressures  and  unburned  flame  gas 
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Fig.  4 SCATTERED  7.IGHT  SICNAL  VS.  BC13  FLOW 


flow  rates,  yielded  qualitatively  similar  results;  condensation  was  readily 
induced  in  flames  0P6  and  P6  but  was  not  observed  in  P3,  0P3,  or  0U4.  The 
small  amounts  of  scattered  light  observed  in  the  jets  of  the  latter  flames, 
which  were  reduced  by  using  the  smaller  nozzle,  are  attributed  to  unwanted 
oarticle  formation  and  ablation  of  condensate  in  the  reservoir  section — not  to 
gas  phase  condensation  anywhere  in  the  test  flame  flow  stream. 

2 . Mass  Spectrometry 

Supersonic  jets  produced  by  the  expansion  of  many  of  the  H2/02/ 
N2  flames  of  Table  3 seeded  with  BC13  were  examined  mass  spectrometrically . 

Upon  the  introduction  of  the  BC13  to  H2-rich  flames,  ions  at  masses  35-38  and 
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37),  HC1  (36  and  38)  and/or  HB02  (43  and  44),  and  B(0H)2  (45  and  44).  Whether 
B02  is  present  is  uncertain;  there  is  not  enough  data  to  quantitatively  define 
I the  cracking  patterns  for  HB02  and  B(OH)2  of  which  m/e  = 43  and  42  may  be  part. 

J In  addition  to  these  species,  02-rich  flames  exhibit  small  concentrations  of 

| species  at  masses  51-55.  Those  may  be  due  to  CIO  (51  and  53)  and  H0C1  (52  and 

54)  although  thermodynamic  considerations  predict  very  small  concentrations 
i (<  10-5  mole  fraction)  of  these  species  in  the  flame  and  throughout  the  nozzle 

expansion. 
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Of  the  four  test  systems  of  Fig.  4,  the  jets  produced  by  the  j-j 

two  cooler  flames,  when  supplied  with  amounts  of  BC13  below  that  required  for 
condensation,  exhibit  a wider  variety  of  boron  combustion  products  and  much 
lower  concentrations  of  HB02  than  do  the  hot  flame  jets.  In  addition  to  the 
species  observed  from  the  hot  flames,  ions  appear  at  masses  70-72,  88,  89,  106, 

113-115,  and  130-132.  An  example  of  such  an  ion  spectrum  is  given  in  Fig.  5. 

This  data  was  taken  at  low  mass  spectrometer  resolving  power  to  enhance  sensi- 
tivity and  reduce  mass  discrimination  effects  in  the  mass  filter  but  even  so, 
most  of  the  ions  listed  are  separately  observable. 

The  species  giving  rise  to  the  observed  ions  have  now  been 
identified  in  nearly  every  case  and  some  cracking  patterns  have  been  discerned. 

The  trimer  of  HB02  has  a complex  cracking  pattern  at  high  (*v  70  V)  ionizing 
energies,  giving  rise  to  ions  throughout  the  groups  about  114  and  131.  Many 
of  these  fragment  ions  decrease  in  intensity  or  disappear  rapidly,  however,  as 
ionizing  energy  is  reduced  and  at  28  V,  the  energy  at  which  the  data  of  Figs.  5 
to  7 were  taken,  fragmentation  appears  to  be  relatively  unimportant.  Mass  70 
is  certainly  due  to  B203.  It  is  not  unambiguously  clear  how  the  ions  at  masses 
72,  88,  89,  and  106  are  related  to  one  another  or,  from  what  neutral  species 
all  of  them  are  formed;  from  examination  of  ion  spectra  as  a function  of  ionizing 
voltage,  it  appears  that  88  and  at  least  some  of  the  signal  at  72  are  indicative 
of  (HB02)2  and  that  masses  89  and  106  arise  from  the  adduct  molecule  HB02»H3B03. 

Additional  support  for  these  contentions  is  derived  from  the  fact  that  the 
intensities  of  masses  115  and  132,  on  the  other  hand,  seem  to  be  third  order 
in  added  boron  at  low  ionizing  energy  and  therefore  appear  to  contain  three  B 
atoms;  mass  106  does  not  yield  unambiguously  to  this  analysis  because  of  its 
relatively  low  signal  level. 

It  appears  quite  clear  from  the  data  of  Fig.  6 that  mass  44  is 
first  order  in  added  boron  and  therefore  contains  but  one  B atom;  similar  plots 
(not  shown)  are  obtained  for  mass  45.  The  relative  abundances  of  these  two 
ions  vary  considerably  from  flame  to  flame,  however,  and  do  not  correlate  at 
all  with  the  isotopic  distribution  for  boron  assuming  they  both  n,rise  from  the 
species  H2B02.  Mass  44  is  much  (2-10  fold)  more  prevalent  in  hot  02-rich 
flames  and  mass  45  becomes  predominant  at  lower  temperatures  in  H2-rich  enviren- 
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ments.  The  behavior  of  mass  45  approximately  parallels  that  of  mass  62  which 
is  attributed  to  H3B03  on  the  basis  of  its  observed  first  order  dependence  on 
added  BC13.  The  ratio  of  the  intensities  of  masses  45  and  62,  however,  is  not 
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Fig.  6 MB/MS  RESPONSE  AS  FUNCTION  OF  TOTA' 
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strongly  dependent  on  ionizing  energy  and  if  H2B02  is  an  ion  source  fragment 
of  H3BO3  the  daughter  ion  is  roughly  10  times  more  abundant  than  the  parent 

over  the  ion  source  energy  range  25-70  V.  It  seems  more  likely  that  the  bulk  ] 

\ 

of  mass  45  is  due  to  an  actual  neutral  species  H2B02  rather  than  exclusively  j 

a daughter  ion  of  H3B03.  I 

] 

J 

In  principle  it  is  possible  to  determine  the  number  of  boron  1 

atoms  in  a given  species  by  examining  the  relative  abundances  of  the  isotopes  j 

within  a group  of  ions.  A molecule  containing  a single  B atom  will  produce  j 

two  ions  at  adjacent  mass  numbers,  with  the  higher  molecular  weight  being 
about  four  times  more  intense  than  the  lower;  two  B atoms  yield  three  adjacent  j 

peaks  with  relative  abundances  4:31:66  in  order  of  increasing  mass;  three  B J 


DISTANCE  FROM  NOZZLE  EXIT,  cm 

Fig.  7 MB/MS  RESPONSE  TO  B-CONTAININO  SPECIES  IN  SUPERSONIC  JET 

Reservoir:  Flame  OP6  = 1500  K , P = 150  Torr,  BC13  = 0.6%; 

Jet:  T = 'v  650  K,  P = 4.5  Torr,  velocity  = 2.3  x 10s  cm  sec-1 
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atoms,  four  peaks  in  the  ratios  6:6:86:372:535.  On  this  basis  the  ion  spectrum 
in  Fig.  5 is  clearly  indicative  of  three  B-atom  species  at  masses  131  ± 1 and 
114  ± 1 although  in  both  these  cases  the  smallest  ion  in  each  set  of  four  peaks 
(masses  129  and  112,  respectively)  is  lost  in  the  background  noise.  It  is 
also  apparent  from  Fig.  5 that  two  different  species  are  responsible  for  the 
peaks  at  88  ± 1,  although  here,  too,  the  least  intense  components  in  the  group 
are  not  visible.  Unfortunately  this  analysis  reveals  nothing  about  the  relative 
contributions  of  the  two  possible  sources  of  mass  45;  both  the  molecule  H2B02 
and  the  daughter  ion  of  H3B03  contain  only  one  B atom. 

The  data  of  Fig.  7 were  obtained  by  traversing  the  nozzle  over 
a 6 cm  distance  in  front  of  the  mass  spectrometer  sampling  cone  while  maintain- 
ing the  reservoir/nozzle  exit  pressure  ratio  required  for  perfect  expansion. 

The  burner-to-nozzle  distance  in  the  reservoir  section  of  the  system  was  held 
constant  to  avoid  interference  by  kinetic  effects  in  the  flame  gases  upstream 
of  the  nozzle  inlet.  The  exit  pressure  conditions  required  for  perfect  expan- 
sion of  the  jet  were  easily  determined  by  monitoring  the  mass  spectrometer  ion 
source  pressure  as  the  nozzle  exit-to-sampling  cone  distance  was  varied  as 
shown  in  Fig.  8.  Sizable  excursions  in  the  mass  flow  into  the  ion  source  are 
observed  under  both  overexpanded  and  underexpanded  conditions  due  to  changes 
in  gas  density  produced  by  the  formation  of  shocks  in  the  jet.  As  conditions 
for  perfect  expansion  are  approached,  the  magnitude  of  the  pressure  fluctua- 
tions decreases  and  the  gas  flow  into  the  mass  spectrometer  becomes  nearly  con- 
stant throughout  the  length  of  the  jet.  In  the  flame  jet  of  Fig.  8,  this  con- 
dition is  closely  approximated  at  an  ambient  (expansion  chamber)  pressure 
Poo  = 5.7  Torr  in  fairly  good  agreement  with  the  calculated  value  of  4.7  Torr 
given  in  Table  3 for  fully  frozen  composition.  It  is  also  interesting  to  note 
that  the  apparent  position  of  the  first  shock  in  the  underexpanded  flow  at 
1.5  Torr  is  located  about  5 exit  diameters  (L/D)  downstream  of  the  nozzle. 

This  compares  quite  well  with  a distance  of  5.9  diameters  predicted  from  the 
correlation 

(L/D)  = 0.67  (Pres/Pco)1/2  (6) 

where  Pres  is  the  combustion  chamber  pressure  (150  Torr).  The  nozzle  expansion 
and  jet  characteristics  of  even  the  smaller  of  our  two  nozzles  appear  therefore 
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DISTANCE  FROM  NOZZLE  EXIT,  cm 


Fig.  8 MASS  SPECTROMETER  ION  SOURCE  PRESSURE  AS  FUNCTION  OF  JET  POSITION 

AND  AMBIENT  PRESSURE 

Reservoir:  Flaiue  OP3,  P = 150  Torr 


to  be  fairly  close  to  those  predicted  on  the  basis  of  ideal  behavior;  undesir- 
able boundary  layer  effects  and  complications  due  to  shock  formation  within  the 
nozzle  are  apparently  unimportant.  This  is  a very  important  observation  since 
no  attempt  was  made  in  this  work  to  experimentally  determine  static  temperatures 
in  the  supersonic  jet  and  instead  the  temperatures  were  computed  on  the  basis 
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of  ideal  gas  dynamic  behavior.  The  use  of  this  assumption  appears  justified 
in  view  of  the  above  and  the  jet  temperature  is  therefore  considered  to  be 
known  as  well  as  the  measured  flame  temperatures  in  the  reservoir  section 
(typically,!  40  K) . 

As  the  total  boron  concentration  corresponding  to  the  threshold 
of  condensation  is  approached  in  the  jets  of  the  cooler  flames,  white,  powdery 
boron  deposits  were  observed  to  collect  in  two  critical  spots  within  the 
apparatus — the  mass  spectrometer  sampling  cone  and  the  nozzle  throat.  Attempts 
to  prevent  this  condensation  on  the  sampling  cone  by  electrically  heating  the 
lower  section  and  base  of  the  cone  were  only  partially  successful;  as  particles 
began  to  form  in  the  jet,  orifice  clogging  followed  almost  immediately  and  the 
experiment  had  to  be  interrupted  for  cleaning.  Deposition  of  solids  in  the 
throat  of  the  nozzle  could  be  tolerated  for  longer  periods  of  time.  The 
immediate  effect  of  this  deposition  is  a decrease  in  effective  throat  diameter 
which  alters  the  expansion  ratio  and  mass  throughput.  This  was  reflected  in 
steady,  small  changes  in  the  jet  pressure  and  was  detectable  as  a decrease  in 
the  mass  spectrometer  ion  source  pressure.  In  order  to  minimize  these  prob- 
lems, particularly  near  the  condensation  threshold  where  they  were  most  severe, 
BC13  flow  rates  were  first  carefully  set  with  the  additive  being  dumped  into 
the  vacuum  system  downstream  of  the  expansion  chamber  and  then  diverted  into 
the  flame  for  periods  of  only  a few  seconds  during  which  specific  ion  current 
measurements  were  made. 

From  the  data  of  Fig.  6,  it  is  apparent  that  at  total  boron 
concentrations  approaching  1%,  species  other  than  those  cited  must  be  produced; 
negative  deviations  from  the  direct  integral  order  dependences  indicated  in 
Fig.  6 are  observed  for  all  species  of  mass  132  or  less.  Rapid,  low  resolution 
scans  of  the  upper  portion  of  the  mass  spectrum  (132  < m/e  < 500)  under  these 
conditions  revealed  peaks  at  the  approximate  masses  158,  185,  202,  230,  246, 
274,  300,  315,  385,  and  455.  Efforts  to  determine  the  dependences  of  these 
species  on  added  boron  proved  fruitless;  the  boron  concentration  range  within 
which  they  were  observable  prior  to  orifice  clogging  was  simply  too  small  for 
meaningful  measurements.  The  peaks  identified  were  observed  at  an  ionizer 
energy  of  25  volts,  h’o  significant  changes  in  the  qualitative  features  of  the 
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mass  spectra  were  detected  at  energies  up  to  70  volts  although  very  little 
effort  was  devoted  to  obtaining  the  high  energy  data. 

At  the  lowest  total  boron  concentration  which  yields  detectable 
quantities  of  masses  with  500  > m/e  > 132,  158  and  230  are  most  intense.  How- 
ever, the  ions  within  this  entire  mass  range  comprise  only  a small  fraction  of 
the  total  actually  present.  Turning  off  the  quadrupole  d.c.  voltage  at  a fixed 
mass  number  setting  allows  all  ions  with  masses  equal  to  or  greater  than  that 
value  to  be  transmitted  to  the  detector.  By  comparing  the  relative  magnitudes 
of  the  signals  for  species  above  132  and  above  500,  one  obtains  a measure  of 
the  proportion  of  high  molecular  weight  species  represented  by  the  ions  listed 
above,  which  are  observable  within  the  presently  accessible  mass  range  of  the 
instrument.  These  tests  reveal  that,  at  the  limit  of  detectability  for  ions 
with  132  > m/e  > 500,  the  boron-containing  species  are  partitioned  as  follows 
between  various  mass  ranges:  40  < m/e  < 140  = 30%;  140  < m/e  < 500  = 10%; 

500  < m/e  = 60%.  At  higher  boron  concentrations,  the  higher  molecular  weight 
species  comprise  an  extremely  rapidly  increasing  proportion  of  the  total. 

D,  DISCUSSION 

At  equilibrium,  nearly  all  of  the  boron  present  in  the  jets  of  any 
of  these  flames  would  be  in  the  form  of  liquid  or  solid  B203.  This  is  illus- 
trated by  the  'Complete  Equilibrium'  data  given  in  Table  5 for  two  of  the  four 
tests  systems — the  jets  of  flame  0P6  (which  exhibits  particle  formation)  and 
flame  0P3  (which  does  not).  The  first  attempt  to  interpret  the  results  of  our 
light  scattering  experiments  was  therefore  based  on  a model  involving  homoge- 
neous nucleation  of  B203.  If  it  is  assumed  that  the  gas  phase  reactions  occur- 
ring in  the  nozzle  are  very  fast  and  that  particle  formation  is  limited  by  the 
rate  of  nucleation  and  condensation,  the  distributions  of  B-containing  products 
in  these  two  jets  assume  the  values  labeled  'Equilibrium,  No  Condensation'  in 
Table  5.  The  resultant  supersaturation  ratios,  S,  for  B203  have  also  been 
computed;  they  are  simply  the  ratio  of  the  'No  Condensation'  values  to  those 
at  'Complete  Equilibrium'  (with  liquid  B203).  It  may  be  seen  from  these  num- 
bers that  the  less  supersaturated  of  the  two  systems  given  is,  in  fact,  the 
one  which  exhibits  particle  formation — the  direct  opposite  of  the  result  antici- 
pated on  the  basis  of  the  simple  analysis.  This  inconsistency  prevails  in  the 
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Table  5 

DISTRIBUTIONS  OF  B-CONTAINING  PRODUCTS  AFTER  NOZZLE  EXPANSION 
FROZEN  AND  WITH  AND  WITHOUT  CONDENSATION 


(Total  B *=  1 mole  X.  Reservoir  P e 100  Torr 


Complete  Equilibrium 


Flame 

0P6 

OP  3 

T(reservolr)  * 

1500 

1700 

T K 

(nozzle  exit) ’ 

650 

780 

Jet  Mole  Fractions 

HBO, 

1.3(-10) 

7 . 5 (-7) 

BO, 

4. 9(-14) 

9.5  <—12 't 

HjBO, 

1. 3 (— 6 > 

2. 7 (-4) 

(HBO,), 

5.0(-8) 

9. 3(-5) 

B,0,(?J 

1.0(-2) 

1.0(-2) 

Supersaturation 

Ratio 

1.0 

1.0 

Equilibrium,  Frozen, 


No  Condensation 

Chamber 

Composition 

OP  6 

OP  3 

0P6 

OP  3 

1500 

1700 

1500 

1700 

650 

780 

650 

780 

7 . 7 (-9) 

3.6(-6) 

1.0C-2) 

1 . 0 (-2 ) 

2 . 9(-12) 

4 . 6 (-11 ) 

9. 2 (-6) 

3.  9(-5) 

7 . 8(-5) 

1.1 (-3) 

2 . 3 (-6) 

5. 9 (-7) 

1.0(-2) 

8.91-3) 

5.7  C-9) 

3.4 (-11) 

3.3 

21 

2.3(18) 

1.1(10) 

instances  of  the  other  two  flame  jets  of  Fig.  4;  corresponding  supersaturation 
values  for  flame  jets  P6  and  0U4  are  8.2  and  27,  respectively.  Thus,  the  pair 
of  systems  with  lowest  S (3.3  and  8.2)  produce  particles  and  the  two  of  highest 
S (21  and  27)  do  not. 


An  alternative  possibility  is  that  the  nozzle  chemistry  is  not  fast, 
the  extreme  case  in  this  instance  being  a nozzle  exit  composition  frozen  at 
the  flame  reservoir  values.  The  mole  fractions  data  in  the  'Frozen  Chamber 
Composition'  columns  of  Table  5 are  indicative  of  the  jet  compositions  in  this 
situation  and  it  may  be  seen  that  the  partial  pressures  of  B203  are  large  com- 
pared to  either  of  the  other  two  cases.  This  abundance  is  reflected  in  the 
extremely  large  values  of  S which  result.  For  the  two  flame  jets  P6  and  0P6, 

S = 6.7  x 1018  and  2.3  x 1018,  respectively;  for  0P3  and  0U4,  S = 1.1  x 10l°  and 
3.0  x 1016.  In  these  frozen  composition  cases  then,  the  relative  super- 
saturation ratios  are  at  least  in  the  same  order  as  the  approximate  tendency 
to  exhibit  condensation,  i.e.  P6  > 0P6  > 0P3  or  0U4.  Unfortunately,  the  S 
values  are  so  large  as  to  be  meaningless  with  respect  to  assessing  critical 
nucleus  size.  If  it  is  assumed  that  S = 1017  represents  the  critical  super- 
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saturation  and  that  the  surface  tension,  G,  of  the  liquid  B203  particle  is 
% 80  dynes  cm-1,  then  according  to  the  Kelvin  equation 


2 a V 


kT  In  S 


(7) 


(where  V is  the  B203  molecular  volume)  the  radius  of  the  critical  nucleus, 
rc  S 2.5  x 10~9  cm  at  700  K.  This  radius  corresponds  to  a critical  nucleus 
volume  of  6.5  x 10  26  ml  or  about  10  3 B203  molecular  volumes,*  indicating 
that  a single  B203  molecule  is  already  larger  than  the  critical  size  for  con- 
tinuous growth. 


Experimentally  (see  below)  the  exit  plane  composition  is  shown  by 
the  mass  spectrometric  results  to  lie  intermediate  between  the  'Frozen'  and 
'Equilibrium,  No  Condensation'  conditions  displayed,  in  Table  5 — at  least  at 
total  B concentrations  where  reproducible  measurements  could  be  made,  viz., 
below  those  which  induce  particle  formation  (but  well  above  saturation) . This 
data  together  with  the  inconsistencies  in  the  results  of  the  analysis  based  on 
simple  homogeneous  nucleation  of  B203  makes  it  impossible  to  predict  the  condi- 
tions prerequisite  to  condensation  from  simple  nucleation  theories.  The  behav- 
ior of  individual  gas  phase  species  as  revealed  by  the  mass  spectrometric 
results  provides  far  better  indications  of  the  mechanism  involved. 


Before  venturing  into  the  nature  of  the  high  molecular  weight  species 
which  may  ultimately  lead  to  particle  formation,  a few  comments  concerning  the 
simpler  B-containing  gaseous  compounds  are  in  order.  The  molecules  HB02,  B203, 
H3B03  and  (HB02)3  are  all  observed  in  abundance.  The  presence  of  none  of  these 
in  the  supersonic  jets  is  surprising;  all  are  well-known  gas  phase  species  for 
which  thermochemical  data  are  available.  The  existence  of  H2B02  (probably 
better  written  as  B(0H)2),  (HB02)2,  H*,B205  and  H2B305  was  not  anticipated, 

however.  Of  these  species  only  B(0H)2  has  been  observed  previously. 14  Unfor- 
tunately, there  is  little  of  a quantitative  nature  that  can  be  done  with  these 
data;  the  systems  of  this  study  in  which  these  new  entities  are  observed  are 


* As  is  customary  in  such  calculations,  the  molecular  volume  (Vra)  is  computed 
from  condensed  phase  macroscopic  density;  i.e.  Vm  = M/pN  where  M = molecular 
weight  (69.6)  and  N is  Avogadro's  number. 
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in  a severe  state  of  disequilibrium  so  it  is  very  perilous  to  attempt  calcula-  j 

tions  of  thermodynamic  properties  and  the  meager  information  on  bulk  gas  free  ] 

radical  concentrations  gives  rise  to  large  uncertainties  in  most  attempted  j 

estimates  of  the  individual  rates  of  formation  or  decay  of  observed  species.  ; 

For  reactions  not  involving  bulk  gas  constituents,  these  uncertainties  dis-  j 

appear,  of  course,  and  it  becomes  possible  to  estimate  rates  of  such  reactions  j 

subject  only  to  the  errors  associated  with  the  quantitative  interpretation  of  j 

the  mass  spectral  data.  J 

The  polymerization  of  HB02  is  unlikely  to  require  the  participation  j 

of  a third  body.  In  all  likelihood,  the  energy  released  in  the  dimerization,  1 

which  is  probably  of  the  order  of  20-30  kcal,  \ 

2HB02  (HB02)2  (8) 

can  be  readily  accommodated  in  the  18  vibrational  and  3 rotational  degrees 
of  freedom  available  for  a time  long  enough  for  collisional  stabilization. 

Molecules  of  comparable  complexity  and  even  more  stringent  energy  dissipation 
requirements  frequently  display  second-order  kinetics  in  this  pressure  range.23 
Therefore  this  reaction  and  the  trimerization 

(HB02)2  + HB02  -+  (HB02)3  (9) 

as  well,  are  probably  not  pressure  dependent*  and  the  data  of  Fig.  7 may  be 
employed  to  deduce  rate  coefficients  for  these  processes  directly  from  the  4 

partial  pressures  of  the  reactants.  A linear  replot  of  the  mass  spectrometer  j 

response  to  these  three  species  as  functions  of  reaction  time  (flow  distance/ 
jet  velocity)  is  given  in  Fig.  9.  The  rate  coefficient  for  Reaction  (9),  k9, 
may  be  obtained  directly  by  measuring  the  initial  slope  of  the  (HB02)3  pro-  * 

file  and  dividing  it  by  the  product  [HB02 ) [ (HB02) 2 ] . Implicit  in  this  proce- 
dure are  the  assumptions  that:  (1)  Reaction  (9)  is  the  only  source  of  the 

trimer,  (2)  at  t = 0,  the  rate  of  the  reverse  of  Reaction  (9)  is  negligible 


* This  point  could  be  readily  checked  by  varying  jet  static  pressure.  Unfor- 
tunately since  the  determination  of  these  rate  coefficients  was  not  the  major 
thrust  of  this  effort,  time  did  not  permit  the  execution  of  such  experiments. 

23.  Benson,  S.W.,  The  Foundations  of  Chemical  Kinetics  (McGraw  Hill,  New  York, 
1960),  Chaps.  XI  and  XII. 
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The  left-hand  side  of  this  expression  may  be  obtained  directly  from  Fig.  9. 
Then  adopting  assumption  (3),  above,  concerning  the  concentrations  and  ion 
current  relationships  and  using  the  computed  local  values  of  kg  one  obtains 
values  for  k8  throughout  the  jet.  The  constancy  of  these  results  over  the 
range  of  flow  times  for  which  data  are  available  is  a qualitative  measure  of 
the  validity  of  the  mechanistic  model.  If  other  processes  intervened  signifi- 
cantly, it  would  be  extremely  fortuitous  if  the  derived  values  for  k8  were 
constant  over  an  appreciable  distance  or  time.  In  fact,  these  calculations 
produce  numbers  with  positive  and  negative  deviations  no  larger  than  15%  of 
the  average  for  the  first  1.2  x 10-5  sec  of  flow  time  after  which  a distinct 
decrease  in  the  results  is  observable.  The  latter  effect  is  probably  due  to 
either  or  both  reactions  approaching  equilibrium  and  the  onset  of  either  or 
both  of  the  reverse  reactions.  The  values  for  these  two  rate  coefficients  at 
the  early  flow  times,  where  they  should  be  most  accurate,  are  given  by 

k8  = (8.1  ± 1.0)  x 10-11  ml  molecule-1  sec-1  (11) 

k9  = (3.0  ± 1.0)  x 10-1°  ml  molecule-1  sec-1  (12) 

at  a jet  temperature  of  650  K and  an  average  [HB02]  % 2 x 10-1'*  ml-1. 


The  absence  of  higher  polymers  of  HB02  from  the  ion  spectrum  under 
any  conditions  seems  to  indicate  that  the  trimer  formation  is  a mechanistic 
dead-end.  This  is  also  inferred  from  the  persistence  of  a third  order 
dependence  of  (HB02)3  on  a'dded  boron  at  relatively  large  addition  rates;  many 
other  species  including  B203  , H3B03  and,  most  notably,  the  entity  with  mass 
115  (H<,B205)  exhibit  deviations  from  integral  value  dependences  on  additive 
at  concentrations  significantly  below  the  points  at  which  (HB02)3  begins  to 
deviate;  this  is  illustrated  by  the  points  corresponding  to  mass  115  in  Fig.  6. 
Behavior  of  this  kind  sets  in  as  chemical  reactions  which  act  as  a sink  for 
the  species  in  question  become  important;  the  kinetic  expression  for  their 
steady-state  concentration  is  then  no  longer  simply  a function  directly  propor- 
tional to  precursor  concentrations.  The  (HB02)2  (88)  curve  in  Fig.  6 illus- 
trates this  point  in  that  it  departs  from  second  order  dependence  on  HB02  at 
additive  levels  where  (HB02)3  formation  becomes  a significant  sink  via 
Reaction  (9),  i.e.,  when  the  second  term  on  the  right-hand  side  of  Eq.  (10) 
becomes  significant. 
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A corollary  of  this  argument  is  that  those  species  which  exhibit 
deviations  from  simple  integral  kinetic  dependences  at  the  lowest  additive  con- 
centrations are  the  species  which  undergo  further  reaction  most  readily.  It 
is  among  these  species  that  the  mechanistic  pathway  to  the  high  molecular 
species,  which  comprise  'pseudo-condensation  nuclei',  must  be  sought.  Further 
clues  as  to  the  identity  of  these  critical  nucleation  precursors  may  be  obtained 
from  the  molecular  weights  of  the  species  observed  at  high  (near-condensation 
onset)  total  boron  concentrations.  Considering  the  abundant  simpler  molecules 
HB02,  H3B03  and  B203  as  monomers  in  the  vocabulary  of  classical  nucleation 
schemes,  one  may  subtract  the  respective  molecular  weights,  44,  62  and  70  from 
the  observed  higher  mass  species  to  deduce  the  identity  of  the  molecule  to 
which  they  adhere.  For  example,  in  the  instance  of  the  polymer  with  mass  158 
to  160  (identified  nominally  as  mass  158  in  Section  III.C.2  above),  the  possible 
origins  are  44  + 115  (HB02  + H2B305),  70  + 88  (B203  *+  (HB02)2),  and  63  + 95. 
Since  the  low  molecular  weight  spectrum  does  not  reveal  any  constituent  with 
mass  95  and  since  the  dimer,  (HB02)2,  appears  to  be  nicely  accounted  for  in 
the  trimer  formation  scheme,  the  first  of  the  three  possibilities  appears  by 
far  the  most  likely.  Reinforcing  this  hypothesis  is  the  observation  that 
H2B305  (115)  begins  to  deviate  from  linear  behavior  in  Fig.  6 relatively  early 
indicating  its  significant  participation  in  some  reaction  which  depletes  its 
concentrations. 

Analogous  arguments  when  applied  to  other  observed  polymeric 
molecules  lead  to  the  mechanistic  hypothesis  partially  outlined  in  Fig.  10. 

This  reaction  scheme  accounts  for  all  the  species  observed  in  this  work  but  is, 
by  no  means,  the  only  possible  route  for  their  formation.  For  example,  at 
some  stage  in  the  polymerization  the  interactions  between  two  entities  each 
containing  more  than  three  B atoms  surely  become  important  and  the  reactions 
postulated  will  be  bypassed.  Similarly,  the  reactions  written  are  not  meant 
to  uniquely  and  exclusively  define  the  elemental  steps  involved;  for  example, 
the  addition  of  HB02  as  written  in  several  instances  may  well  be  the  combined 
result  of  H3B03  addition  accompanied  or  followed  immediately  by  elimination  of 
a molecule  of  H20.  The  intermediate  adduct  in  this  case  would  be  too  unstable 
to  survive  electron  impact  ionization  or  be  of  such  short-lived  duration  as  to 
render  its  steady  state  concentration  too  low  for  detection. 
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Fig.  10  PARTIAL  HYPOTHETICAL  POLYMERIZATION  MECHANISM 


However,  despite  these  uncertainties,  the  picture  which  emerges  from 
this  analysis,  strongly  indicates  that  there  is  no  distinct  step  that  can  be 
identified  with  dehydration  and  that  no  single  monomeric  species  can  be  iso- 
lated up  m which  to  base  a simple  model.  As  the  scheme  in  Fig.  10  progresses 
from  left  to  right,  the  proportion  of  boron  in  the  observed  species  gradually 
increases  implying  the  progressive  elimination  of  H20.  The  fact  that  a series 
of  polymers  separated  by  18  mass  units  is  not  observed  would  seem  to  indicate 
that  the  dehydration  occurs  predominant Jy  among  the  lower  molecular  weight 
species,  i.e. , via  the  simple  reactions  interconnecting  HB02 , H3B03,  B203,  etc. 
It  is  possible,  though,  that  the  members  of  these  series  other  than  those 
observed  are  simply  present  at  concentrations  too  low  for  detection.  The  model 
is  unfortunately  not  yet  complete. 
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IV,  CONCLUSIONS  AND  RECOMMENDATIONS 

It  has  been  demonstrated  that  the  energetically  important  condensation 
of  boron  combustion  products  cannot  be  accurately  predicted  from  models  based 
on  simple  homogeneous  nucleation  theory.  Mass  spectrometric  examination  of 
systems  actually  undergoing  nucleation  reveal  that  the  reaction  mechanisms  by  I 

which  HB02  is  converted  to  condensed  phase  B203  are  extremely  complex  and  that 
a number  of  'monomeric' species  are  involved.  The  essential  process — that  of  ! 

the  elimination  of  H20 — occurs  gradually  rather  than  at  any  discreet  point  or 
) ; over  any  limited  range  of  polymer  molecular  weights.  A complex  array  of  inter- 

mediate polymeric  species  with  widely  variable  H/O/B  proportions  is  present  in 
these  dynamic  condensing  combustion  product  gas  streams.  The  reactions  leading  ) 

to  the  formation  of  the  dimer  and  trimer  of  HB02,  although  readily  anticipated 
on  intuitive  and  thermochemical  grounds  do  not  appear  to  be  essential  steps  in 
the  overall  process  leading  to  the  production  of  condensed  phases.  Rate  coeffi- 
cients for  these  two  reactions,  viz., 

2HB02  (HB02)2 

and  HB02  + (HB02)2  + (HB02)3 

have  been  measured  here  at  a temperature  of  650  K as  8 x 10-11  and  3 * 10-1° 

ml  molecule-1  sec-1,  respectively.  However  the  portion  of  the  overall  mecha- 

* 

nism  which  they  comprise  appears  to  be  incidental  in  the  overall  scheme;  as 
nucleation  proceeds,  these  dimers  and  trimers  may  well  become  involved  but 
other  reactions  appear  to  play  a much  more  significant  role  in  the  production 
of  higher  molecular  weight  species. 

It  is  quite  likely  that  complicated  phenomena  of  the  sort  encountered  in 
the  study  are  operative  in  a great  many  other  practical  systems  as  well.  There 
are,  in  fact,  very  few  real-life  situations  to  which  one  may  apply  conventional 
nucleation  theories  with  much  prospect  of  attaining  accurate  predictions  due 
to  an  almost  certain  knowledge  of  the  complexity  of  the  chemistry  involved. 

The  work  performed  on  this  contract,  although  rt:  does  not  yet  provide  a com- 
plete picture  of  boron  oxide,  nucleation,  demonstrates  that  the  chemistry  of 
such  processes  can  be  explored  in  great  detail  even  in  situations  where  the 
overall  chemical  composition  is  quite  complex. 
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To  complete  the  effort  initiated  with  this  study,  the  experimental 
capability  of  the  molecular  beam  mass  spectrometer  should  be  extended  to  higher 
mass  numbers  to  facilitate  examination  of  the  later  stages  of  nucleation  and 
particle  growth.  This  work  should  also  include  an  investigation  of  the  effects 
of  additives  or  potential  fuel  impurities  which  might  accelerate  the  nuclea- 
tion and/or  particle  growth  processes.  Candidate  'catalysts'  include  charged 
species  which  may  induce  heterogeneous  ionic  nucleation  and/or  more  easily 
condensed  materials  such  as  silicates  which  may  serve  as  nuclei  for  boron 
combustion  products  by  virtue  of  their  mutual  liquid  phase  solubilities.  If 
H any  of  these  effects  prove  significant  it  may  be  possible,  through  the  inclu- 

j sion  or  elimination  of  a relatively  small  amount  of  added  material  to  effec- 

tively control  particle  and  smoke  formation  via  alteration  of  condensation 
kinetics.  The  instrumental  aspects  of  this  extended  work  are  straightforward; 
however,  several  system  modifications  including  jet’ expansion  to  lower  pres- 
I sures  and  enlargement  of  the  first  stage  sampling  capacity  are  indicated. 


The  systems  amenable  to  investigation  via  the  techniques  employed  in  this 
work  are  by  no  means  limited  to  supersonic  jets  containing  metal  oxides.  In 
principle,  a wide  variety  of  other  nucleation  problems  of  considerable  practical 
importance  could  be  attacked  including  those  associated  with  the  formation  of 
atmospheric  aerosols,  condensation  induced  by  the  interaction  of  more  conven- 
tional exhaust  products  with  air  and  smoke  formation — whether  intentional  or 
accidental.  The  essential  experimental  features  of  the  approach  to  be  taken 
in  all  of  these  cases  would  be  the  same  as  those  employed  in  the  present  work. 
Specific  sampling  system  parameters,  mass  spectrometer  response  characteristics 
and  the  design  of  the  apparatus  in  which  the  nucleation/condensation  processes 
are  induced  would  be  tailored  to  the  particular  needs  of  the  study  to  be 
undertaken. 
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The  rate  coefficient  for  the  dissociative  electron  attachment  e + HBOj  -*  BOj  + II  has  been  measured  in  a number  of 
well  characterized  laboratory  test  flames  at  total  pressures  of  1.0  and  0.1 3 atm.  The  data  are  consistent  with  the  Arrhenius 
expression  k = 3.2  X 10-10  exp(-l  1000/71  ml  molecule-1  s-1  between  1730  and  2250  K.  The  failure  of  the  observed 
BOj  formation  rates  to  correlate  in  any  systematic  way  with  hypothetical  three-body  rate  coefficients  demonstrates  that 
three-body  attachment  processes  are  negligible  in  these  flames  and  that  dissociative  attachment  dominates. 


1.  Introduction 

In  recent  years  a number  of  thermochemical 
studies  have  shown  that  the  oxides  and  oxyiiydroxides 
of  boron  [1,2]  and  a number  of  transition  metals 
[3-5]  have  electron  affinities  as  great  or  greater  than 
those  of  the  halogens.  The  abilities  of  these  metal- 
containing  compounds  to  form  stable  negative  ions  at 
high  temperatures  make  them  of  considerable  interest 
in  combustion  plasmas  such  as  MUD  channel  flows 
and  rocket  exhausts  where  the  concentrations  of  free 
electrons  are  important.  Although  these  studies  have 
resulted  in  the  generation  of  the  thcrmochcmical 
data  necessary  to  assess  the  equilibrium  effects  of 
these  electrophilic  compounds,  kinetic  information 
is  still  lacking.  The  work  described  here  comprises 
the  first  measurements  of  the  rate  coefficient  for  the 
reaction 

e + HB02  B02  + H . (1) 

The  measurements  were  performed  in  a number  of 
well-characterized  laboratory  combustion  flames 
spanning  the  temperature  range  1730-2250  K at 
pressures  of  1 .0  and  0. 1 3 atm  ( 1 00  torr). 

2.  Experimental 

Due  to  the  possibility  of  interference  in  the  meas- 


urements of  A:  j by  three-body  electron  attachment 
reactions  it  was  deemed  desirable  to  perform  meas- 
urements of  BOJ  formation  rates  in  both  1 atm  and 
low  pressure  flames.  The  1 atm  flames  used  are  the 
same  H2/02/N 2 systems  employed  in  prior  thermo- 
dynamic work  [1-6].  The  concentration  profiles 
for  all  the  important  neutral  constituents  (I I20,  H2, 

H,  Oil)  as  functions  of  flow  or  reaction  time  in  these 
flames  are  thus  available  from  these  and  other  previous 
studies  and  no  further  work  of  this  nature  was  requir- 
ed at  1 atm.  Typically,  the  results  of  these  concentra- 
tion measurements  span  reaction  times  from  0.5  to 
5 ms  - an  interval  adequate  for  the  boron-containing 
species  of  this  study  to  reach  near-equilibrium  con- 
centrations. Unfortunately,  very  little  neutral  species 
concentration  data  has  been  generated  for  H2-rich 
low  pressure  flames  and  it  was  necessary  to  charac- 
terize these  flames  with  respect  to  their  H-atom  pro- 
files and  temperatures. 

The  low  pressure  measurements  were  performed  on 
an  annular  burner  (inner  flame  diameter  = 2.6  cm, 
overall  flame  diameter  = 4.5  cm)  contained  in  a vacu- 
um housing;  both  burner  sections  are  supplied  with 
H2/02/N2  mixtures  at  flow  rates  of  ^30  ml(STP) 
s-*  per  cm2  of  burner  area.  The  resulting  flames  are 
stable  in  the  pressure  range  50- 1 50  torr.  As  in  prev- 
ious work,  potassium  salts  are  supplied  to  the  inner 
flame  with  an  aspirator  in  the  N2  line.  The  aspirator 
is  operated  at  or  near  1 atom,  which  pressure  is 
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maintained  by  a critical  flow  orifice  located  in  the 
burner  feed  line.  The  boron  additive  compounds  arc 
added  either  from  a hypodermic  infusion  pump  or 
from  the  aspirator  as  a co-solution  with  the  potassium 
salt. 

A convenient  operating  pressure  for  these  flames 
was  found  to  be  about  100  torr.  Such  a pressure  is 
high  enough  to  make  diffusional  ion  transport  unim- 
portant* relative  to  convection  and  low  enough  to 
provide  the  spatial  resolution  necessary  to  obtain 
precise  concentration  profiles.  With  0.1  molar  solu- 
tions of  KNO3  and  H3BO3  in  the  aspirator,  flames 
are  obtained  at  100  torr  containing  total  potassium 
and  boron  concentrations  ( [K|c  and  |B|C,  respec- 
tively) ^ 5 X 101 1 ml" 1 . A small  (<  IT.  but  unmeas- 
ured) amount  of  hydrocarbon  is  added  to  accelerate 
ion  formation.  Electron  concentration  measurements 
are  performed  using  a large  microwave  resonant  cavity 
(7]  especially  constructed  to  accommodate  the  low 
pressure  burner  system.  Cavity  data,  taken  in  the 
down-stream  regions  of  the  flame  where  axial  electron 
concentration  gradients  are  small,  serve  as  a calibra- 
tion for  electrostatic  probe  determinations  [8[  which, 
in  turn,  provide  highly  spatially  resolved  measurements 
of  overall  charged  species  concentrations. 

Ion  identification  and  changes  in  relative  ion  con- 
centrations with  reaction  time  are  monitored  mass 
spectrometrically.  The  ion  mass  spectrometer  is  used 
in  the  same  configuration  [3]  as  in  previous  studies. 

* Based  on  observed  profiles  for  BOi  and  Cl"  in  a 1 00  torr 

flame,  axial  diffusion  velocities  for  negative  ions  are  about 
1 cm  s"1  compared  to  a convective  flow  velocity  of  1500 
cm  s'  ; radial  profiles  have  not  been  taken  but  in  the  first 
8 cm  or  so  of  flow  where  attachment  occurs,  radial  diffu- 
sion may  quite  safely  be  ignored. 

Table  1 

Low  pressure  (100  torr)  flame  characteristics 


The  only  modification  necessary  is  enlargement  of 
the  sampling  orifice  size  when  100  torr  flames  are 
sampled.  At  1 atm  the  sampling  cone  is  pierced  by  a 
0.05  mm  diameter  hole;  at  100  torr,  a cone  with  a 
0.125  mm  diameter  orifice  is  used.  The  axial  burner- 
to-sampling  orifice  distance  is  continuously  variable 
from  0 to  25  cm. 

Flame  temperatures  at  100  torr  are  measured  via 
the  conventional  Na-D  line  reversal  method.  The 
results  for  the  flames  selected  for  use  in  this  study  are 
given  in  table  1 together  with  H-atom  data  at  a num- 
ber of  positions  in  the  flame.  Temperatures  were 
measured  only  in  a region  6-  10  cm  downstream  of 
the  burner  to  avoid  errors  arising  from  chemi-excita- 
tion  or  nonequilibrium  energy  transfer  processes  in 
the  flame  front. 

The  H-atom  concentrations  in  table  1 were  derived 
from  mass  spectrometric  measurements  of  the  relative 
concentralions  of  IIM0O4  and  M0O3  produced  upon 
the  addition  of  small  quantities  (<  10~8  mole  frac- 
tion) of  molybdenum  salts  from  the  aspirator  feed 
system.  These  two  ionic  species  are  linked  through 
the  reaction  (4J 

H + HM0O4  =Mo03  +H20  (2) 

for  which  the  equilibrium  constant  is  given  by  the 
expression  A' 3 - 0.85  exp(  J 0400/7").  Ion-molecule 
reactions  of  this  ivpe  are  sufficiently  rapid  to  be  in 
equilibrium  at  all  points  in  the  test  flames  and  thus 
the  ion  current  ratios  are  a good  indication  of  the 
local  [H]/[H-)0|  ratio.  The  values  for  this  ratio  are 
then  used  together  with  the  equilibrium**  HsO  con- 

**  All  thermodynamic  data  for  equilibrium  calculations  arc 

taken  from  ref.  [9|. 


h2/o2/Nj 

Velocity 

|ll)  (10 

15  ml"') 

unburned  mole  ratios 

T(K) 

(m  s'* ) 

distance  from  burner  (cm) 

2 

4 

6 

8 

10 

12 

14 

25:1.0:2.0 

2100 

11.6 

12.5 

12.5 

12.5 

12.0 

11.4 

11.0 

10.8 

2. 5:1. 0:3.0 

1980 

11.2 

11.5 

11.0 

10.8 

10.5 

10.3 

10.0 

9.6 

3.4: 1.0: 2.0 

1830 

12.1 

12.0 

11.4 

10.7 

10.3 

9.8 

9.4 

90 

3.5: 1.0:3. 0 

1730 

11.7 

11.1 

10.1 

9.6 

9.2 

8 8 

8.5 

8.2 

I 
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centration  to  compute  [H) . Even  though  ronequiii- 
brium  effects  lower  the  actual  [H20]  to  slightly  be- 
low the  equilibrium  value,  the  error  involved  is  small 
since  H20  is  present  in  such  relatively  large  quantities. 

As  a check  on  the  validity  of  the  H-atom  profiles, 
optical  absorption  spectrometry  was  used  to  measure 
[Oil],  The  reaction 

H2+0H-H  + H20  (3) 

is  known  to  he  very  close  to  equilibrium  in  flames  of 
this  sort  and  thus  a profile  of  relative  [H]  can  he  con- 
structed from  measurements  of  relative  [Ollj  and 
well-established  equilibrium  constant  for  reaction  (3). 
No  attempt  was  made  in  these  determinations  to 
establish  absolute  values  for  [Hj;  the  measurements 
instead  serve  as  a check  on  the  mass  spectrometric 
method  above.  In  all  cases,  the  H-atom  decay  curves 
arrived  at  via  both  methods  are  in  quite  good  agree- 
ment and  the  mass  spectrometrically  derived  numbers 
were  adopted  for  use  in  the  kinetic  calculations. 

3.  Results  and  discussion 

Examples  of  the  data  used  for  the  derivation  of  the 
rate  constant  of  reaction  ( 1 ) are  given  in  figs.  I and  2 
taken  at  100  and  760  torr,  respectively.  It  may  be 
noted  that  in  both  cases  the  electron  and  H-atom 


Fig.  1.  Concentration  profiles  in  1980  K (IIj/Oj/Nj  - 
2.S/1. 0/3.0)  100  tort  flame.  |HB02|  = 7 X 1010  raP1 ; 
velocity  = 1 1 .25  m s"1 . 


Fig.  2.  Concentration  profiles  in  2250  K (H2/O2/N2  = 

3.5/1. 0/3.0)  1 atm  flame.  (HB021  = 3 X 10"  ml'1;  velocity 
= 41.6  m s*1 . 

concentration  curves  arc  given  in  absolute  units  but 
that  the  [B02  ] curve  is  expressed  in  realtive  terms 
as  ion  current.  This  is  due  to  the  difficulties  associated 
with  absolute  calibration  of  the  mass  spectrometer. 
Therefore,  in  applying  the  data  to  calculations  of  the 
rate  coefficient,  the  mass  spectrometer  sensitivity 
appears  as  an  additional  unknown. 

It  is  also  necessary  in  deriving  an  expression  for 
the  attachment  rate  constant  to  include  the  effects 
of  the  reverse  of  reaction  ( 1 ).  Assuming  that  the  only 
source  of  B02  is  the  direct  attachment 

e + HB02  -*  B02  + H (1) 

yields  the  following  expression  for  the  formation  of 
B02  : 

d[B02]/df  = k,[e)[HB02l-A-.  JBOJJfHJ  . (4) 

By  further  assuming  that  A'f  = k}/k_ ( , where  AT j 
[=  1 500  cxp(  10000/73]  is  the  equilibrium  constant 
for  the  reaction  and  inserting  / S for  [B02  ] where 
/ is  the  B02  ion  current  (amperes)  and  S is  the  mass 
spectrometer  sensitivity  (molecules  ml”1  ampere”1 ) 
one  may  derive  the  following  equation  of  k j : 

(5, 
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Now,  by  inserting  experimental  numerical  values  of 
[H],  [e],  [HB02]  = [B]c  and  d(/'S)/dr  from  plots 
such  as  figs.  1 and  2 at  various  points  along  the  flow 
stream,  a number  of  simultaneous  equations  for  k j 
and  (/~S)  are  derived.  Pairs  of  these  equations  are 
then  solved  for  the  two  unknowns  and  the  results 
examined  for  consistency.  If  the  calculations  yield 
constant  values  for  both  variables  over  the  entire 
range  of  measurement,  it  is  assumed  that  the  mass 
spectrometer  sampling  system  through-put  did  not 
change  during  the  course  of  the  experiment  and  that 
the  rate  constant  data  is  therefore  valid;  it  is  then 
averaged  with  other  results  for  . The  rate  constants 
obtained  from  all  seven  test  flames  are  given  in  table  2 
together  with  the  flame  temperatures  and  pressures. 
Each  of  the  k\  values  represents  the  average  of  several 
(3-5)  determinations  at  varying  concentrations  (irole 
fractions  between  10“ 9 and  10~6)  of  alkali  metal 
and  total  boron. 

A plot  of  the  bimolecular  rate  constant  data  (kl ) 
versus  7~~’  is  given  in  fig.  3.  The  straight  line  drawn 
through  the  points  is  a least-squares  line,  the  slope  of 
which  yields  an  activation  energy  of  about  92  kJ.  The 
Arrhenius  expression  for  then  becomes 

kl  = 3.2  X 10'10exp(-11000/r) 

ml  molecule-1  s-1  . (6) 

Since  reaction  (1)  is  endothermic  by  78  kJ  at  0 K, 
the  value  92  kJ  appears  to  be  a quite  reasonable  acti- 
vation energy. 

The  possible  importance  of  three-body  attachment 
via  the  reaction 


Table  2 

Average  rate  coefficients  for  e + HBO2  -*  BOj  + H 


TTtorr) 

H2/O2/N2 

unburned 

7TK) 

*12  (ml  molecule 

100 

3. 5/1. 0/3.0 

1730 

8.6<-13)») 

100 

3.4/1. 0/2.0 

1830 

7.0(-  13) 

100 

2.5/1. 0/3.0 

1980 

l.K-12) 

100 

2. 5/1. 0/2.0 

2100 

2 . 8 ( - 1 2 ) 

760 

3.5/1. 0/6.0 

1815 

1.0f-12) 

760 

4. 0/1. 0/4.0 

2020 

t.2(-12) 

760 

3. 5/1. 0/3.0 

2250 

2.0(- 1 2) 

*)  8.6( 

13)  = 8.6  X 10'13 

□ 100  To  rr 
O 760  Torr 


Fig.  3.  Arrhenius  plot  for  c + HBO2  -*  BOj  + H. 

e + B02  + M - B02  + M (7) 

can  be  assessed  by  comparing  the  relative  magnitudes 
of  the  ion  formation  rates  in  the  early  stages  of  the 
flames  where  the  reverse  of  reactions  (1)  and  (7)  are 
least  important  and  the  negative  terms  in  the  kinetic 
expressions  can  be  ignored.  If  the  three-bodv  process 
is  to  have  a significant  effect  on  the  overall  negative 
ion  formation  process,  it  must  be  comparable  in  mag- 
nitude to  the  dissociative  attachment  rate;  thus 

*,  [e][HB02]  <*7[B02][e)[M]  • (8) 

By  miking  the  reasonable  assumption  that  the  reac- 
tion 

H + HB02  - H2  + B02  (9) 

is  balanced  (i.e.,  is  in  equilibrium)  throughout  the 
burned  gases  where  the  attachment  rates  were  deter- 
mined, eq.  (8)  can  be  rewritten 

*,[H2]/K9[H][M]  <*7  , (10) 

where  Af9  is  the  equilibrium  constant  for  reaction  (9). 
Limiting  values  of  the  three-body  attachment  coeffi- 
cient, k-j,  calculated  for  eq.  (10)  arc  given  in  table  3. 

It  can  be  seen  from  this  data  that  if  reaction  (7)  is  to 
be  important  in  the  flames  of  this  study,  its  rate  coeffi- 
cient must  be  very  large  - a factor  of  about  103 
greater,  for  example,  than  either  the  measured  values 
[10,11]  for 
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Table  3 

Lower  limiting  values  for  trial  three-body  rate  contants  (cm6  s~l) 


Hj/Oi/Nj 

unburned 

T(K) 

*7 

k/T 

k/T15 

k/T 2 

3.5/1. 0/3.0 

1730 

6.8(— 27) 

3.688{  — 30) 

8.866(-32) 

21326-33) 

3.4/1. 0/2.0 

1830 

4.3(— 27) 

2.374(-30) 

5.550(— 32) 

1 .2976—33) 

2.5/1. 0/3.0 

1980 

1.9(-27) 

9.404(-  31) 

2.1 13(-32) 

4.7  50(  — 34 ) 

2.5/1. 0/2.0 

2100 

4.2(— 27) 

2.016C-30) 

4.400(-30) 

9600(-34) 

3.5/1. 0/6.0 

1815 

6.4(-27) 

3.50K-30) 

8.2 19(— 32) 

1.929(33) 

4.0/1. 0/4.0 

2020 

2.7(— 27) 

1.309(-30) 

2.9 1 3( —32) 

6.482(-34) 

3.5/1. 0/3.0 

2250 

1 .26—27) 

5.938C— 3 1 ) 

1.2526-32) 

2.639(  34) 

e + 02  + 02  02  + 02  (10 

or  the  calculated  [12]  upper  limit  for  this  reaction. 
The  few  determinations  of  three-bodv  attachment 
rates  currently  available  [13]  typically  lie  in  the 
10'29-10-3°  cm6  s~*  range  with  an  inverse  tem- 
perature dependence  - again  at  variance  with  any  of 
the  three-body  rate  constants  in  table  3.  Moreover, 
the  data  of  table  3 exhibits  scatter  over  nearly  a fac- 
tor of  ten  and  no  systematic  variation  with  tempera- 
ture is  apparent  over  the  range  7fl  to  T~-  It  is  there- 
fore considered  quite  safe  to  assume  that  the  operative 
electron  attachment  mechanism  in  these  flames  is  the 
dissociative  attachment,  reaction  ( 1 ),  and  that  the 
three-body  process,  reaction  (7),  is  unimportant. 
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ABSTRACT 


This  work  comprises  the  first  direct  determination  of  the  kinetics  of  Cl- 
formation  in  flames.  The  experiments  involve  mass  spectrometric  measurements 
of  negative  ion  concentrations  as  functions  of  flow  time  in  the  burned  gases 
of  well-characterized  laboratory  test  flames  seeded  with  small  quantities  of 
potassium  and  chlorine  compounds.  The  results  are  consistent  with  the  two-body 
dissociative  attachment  mechanism 

e + HC1  Cl"  + H 

with  a rate  coefficient  well  represented  by  the  Arrhenius  expression 


k = 3 x 10" lo  exp(-13000/T)  ml  molecule-1  sec-1 

over  the  temperature  range  1730  K to  2475  K.  The  flames  are  operated  at  pres- 
sures of  100  and  760  Torr.  Electron  concentrations  were  measured  using  a micro- 
wave  resonance  cavity  method  and  positive  ion  profiles  (total  charged  species) 
were  determined  with  electrostatic  probes.  H-atom  concentration  profiles 
required  to  take  account  of  the  reverse  detachment  process  are  well  known  for 
the  760  Torr  flames;  in  the  low  pressure  flames  they  were  obtained  using  a mass 
spectrometric  negative  ion  tracer  technique  involving  the  known  ion-molecule 
chemistry  of  molybdenum  additives  and  confirmed  via  optical  absorption  measure- 
ments of  [OH].  Contributions  to  Cl-  formation  from  three-body  attachment  pro- 
cesses are  not  observed.  The  measured  values  of  k are  found  to  be  somewhat 
lower  than,  but  in  general  agreement  with,  those  inferred  from  electron  beam 
and  swarm  studies  but  are  much  lower  than  the  values  previously  derived  from 
measurements  of  the  reverse  reaction  rate. 
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INTRODUCTION 


Although  the  reaction 

e + HC1  -*■  Cl~  + H (1) 

represents  the  most  important  sink  for  electrons  in  a great  many  combustion  sys- 
tems and  plasma  devices,  no  direct  measurements  of  its  rate  coefficient  have 
yet  been  reported.  Electron  beam1’2  and  swarm3  technique  cross  section  measure- 
ments indicate  that  the  rate  constant  is  approximately  given  by  ki  = 10-1° 
exp(l  x 10*/T)  ml  molecule-1  sec-1.  The  rate  coefficient  of  the  reverse  reac- 
tion has  been  measured*  near  room  temperature  in  a flowing  afterglow  to  be 
1.0  x 10  9 ml  molecule-1  sec-1.  Another  determination, 5 made  by  rapidly 
expanding  flame  gases  heavily  seeded  with  chlorine  and  mass  spectrometrically 
monitoring  the  decrease  in  t Cl  1 as  the  expanding  gases  cool,  gives  a rate 
coefficient  for  the  reverse  reaction  similar  to  that  obtained  in  the  flowing 
afterglow.  If  the  temperature  dependence  of  the  rate  coefficient  of  the 
reverse  of  Reaction  1 (k_t)  is  assumed  to  be  small,  one  is  led  to  expect,  by 
taking  the  product  of  k_i  and  the  equilibrium  constant  for  the  reaction,  that, 
at  2000  K,  for  example,  k!  = 5 x 10-9  ml  molecule-1  sec-1.  This  vaiue,  in 
turn,  assuming  the  endothermicity  of  the  reaction  to  be  equal  to  its  activa- 
tion energy,  leads  to  an  Arrhenius  expression  k = 5 x 10-7  exp(l  x 10*/T)  ml 
molecule-1  sec-1.  This  > 10-fold  discrepancy  in  the  pre-exponential  factor 
presents  a considerable  problem  to  those  who  may  wish  to  include  this  reaction 
in  kinetic  models  of  high  temperature  systems  containing  electrons,  chlorine 
and  hydrogen. 

In  an  attempt  to  resolve  this  ambiguity,  direct  measurements  of  the  rate 
of  Reaction  (1)  were  undertaken  in  a series  of  fuel-rich,  premixed  H2/02/N2 
flames  spanning  a temperature  range  cf  1730  to  2475  K.  Flames  burning  at 


760  and  100  Torr  seeded  with  trace  amounts  of  K and  Cl  were  utilized.  These 
flames  provide  a rapidly  moving  stream  of  hot  gases  containing  HC1  in  which 
electron  concentrations  rapidly  increase  as  K is  thermally  ionized.  From  the 
determination  of  the  difference  between  [Cl-]  and  its  equilibrium  concentra- 
tion the  rate  coefficient  is  obtained.  The  results  yield  values  for  ki  similar 
to  those  inferred  from  the  work  of  Refs.  1-3. 

EXPERIMENTAL 

The  flames  used  in  these  experiments  are  produced  on  annular  Meker-type 
burners  which  provide  an  inner  flame  to  which  additives  may  be  supplied,  sur- 
rounded by  an  outer  identical  sheath  flame.  The  760  Torr  flames  are  the  same 
as  those  employed  in  the  first  study  of  this  series6;  they  have  also  been 
extensively  used  in  thermochemical  investigations  and  their  detailed  character- 
istics are  now  well  known.7-9  The  100  Torr  flames  have  an  inner  flame  diameter 
of  2.6  cm  and  an  overall  diameter  of  4.5  cm.  The  low  pressure  burner  is  con- 
tained in  a vacuum  housing  and  is  supplied  with  H2/02/N2  mixtures  at  rates  of 
about  30  ml(STP)sec-1  per  cm2  of  burner  area.  Chlorine  and  potassium  are 
supplied  to  both  types  of  flames  by  aspirating  solutions  of  KC1  or  NH<,C1/KN03 
mixtures  into  the  N2  feed  line  for  the  inner  flame. 


The  temperatures  of  the  760  Torr  flames  have  already  been  well  established; 
those  of  the  100  Torr  flames  were  measured  using  the  conventional  Na-D  line 
reversal  method.  Electron  concentrations  were  measured  using  microwave 
cavities10  which  were,  in  turn,  used  to  calibrate  electrostatic  probe  measure- 
ments of  the  positive  ion  concentration.11  The  probe  measurements  provide 
highly  spatially  resolved  measurements  of  the  positive  ion  concentration  (and 

i 

hence  of  the  electron  concentration  since  negative  ion  concentrations  arc  in  j 

all  cases  relatively  negligible  in  these  experiments).  H-atom  profiles  are 
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veil  known  for  the  760  Torr  flames;  for  the  100  Torr  flames,  they  were  deter- 
mined in  two  ways.  First  a mass  spectrometric  measurement  of  the  ratio  of 
[HMoO<,-]  and  [Mo03-]  was  made  along  the  axis  of  the  flames  seeded  with  traces 
(i  10-8  mole  fraction)  of  molybdenum  salt.  These  ionic  species  are  linked 
through  the  reaction8 

H + HMoO<T  ? Mo03~  + H20  (2) 

which  is  known  to  be  balanced  in  the  760  Torr  flames.  Secondly,  optical 
absorption  measurements  were  made  to  determine  profiles  of  [OH],  This  concen- 
tration can  then  be  related  to  relative  [H]  by  noting  that  the  reaction 

OH  + H2  - H20  + H (3) 

is  very  nearly  balanced  in  flames  of  this  type.  Since  only  relative  [H]  pro- 
files were  obtained  in  the  latter  measurements,  these  serve  only  as  a check 
on  the  mass  spectrometric  absolute  concentration  profile  measurements.  In  all 
cases  the  shapes  of  the  profile  curves  obtained  using  the  two  methods  were 
similar  and  the  mass  spectrometrically  obtained  values  of  [H]  were  used  in  the 
kinetic  calculations. 

The  electron,  positive  ion  or  negative  ion  axial  profile  measure  lents  were 
obtained  by  raising  and  lowering  the  burner  relative  to  the  appropriate  measur^ 
ing  device;  in  the  case  of  the  760  Torr  flames  the  flames  were  examined 
between  0.2  cm  and  5 cm  from  the  burner.  In  the  case  of  tb  100  Torr  flames, 
scans  were  made  over  a range  of  2 to  14  cm  above  the  burner.  In  Table  I the 
principal  flame  parameters  along  with  the  H-atom  profiles  measured  here  and 
in  earlier  work  are  given. 


RESULTS  AND  DISCUSSION 


If  It  is  assumed  that  the  variation  of  [Cl-]  is  controlled  solely  by 
Reaction  (1)  and  its  reverse,  i.e.,  that  3-body  attachment  is  not  important 
and  that  radial  and  axial  diffusion  are  not  important,  then  the  rate  of  change 
in  [Cl-]  along  the  axis  of  the  flame  is  simply  given  by 


^ [Cl-]  = ki [e] [HC1]  - k-i [Cl-] [H]  (4) 

Assuming  that  Ki  = ki/k-i  where  Ki  is  the  equilibrium  constant  for  Reaction  (1) 
and  letting  S be  the  mass  spectrometer  sensitivity,  then  Eq.  (4)  may  be  written 
in  terms  of  the  mass  spectrometer  Cl-  ion  current,  I-,  and  rearranged  to  give 

kl  = vf  £ (!-s)  v[e][HCl]  - "l  j (5) 


where  Vf  is  the  velocity  of  the  burned  gases  and  z is  the  height  above  the 
burner. 

Equation  (5)  may  be  solved  for  ki  directly  by  taking  the  slopes  of  the 
measured  I'profiles  along  with  values  of  [ e]  and  [ H]  at  pairs  of  points  along 
the  I-  profiles  and  then  eliminating  S from  the  pairs  of  equations.  This 
method  has  proved  to  yield  reproducible  (within  a factor  of  about  3)  results 
from  data  taken  over  a period  of  several  weeks  while  obviating  the  necessity 
of  making  an  absolute  calibration  of  the  mass  spectrometer.  This  last  is  an 
important  consideration,  since,  because  of  kinetic  effects  which  may  be  pres- 
ent in  any  'calibration  flame'  in  which  a source  of  negative  ions  is  present, 

(i)  there  is  no  straightforward  method  of  calibrating  the  flame  ion  mass 
spectrometer  for  negative  ions  to  much  better  than  a factor  of  2 and 

(ii)  from  an  experimental  point  of  view,  the  calibration  factor  is  different 
for  each  flame  and  subject  to  drift  from  day  to  day. 
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From  each  I"  curve  three  pairs  of  points  were  typically  selected  and  values 
for  ki  obtained.  In  most  cases  these  values  for  ki  were  in  good  agreement  with 
one  another.  Since  slight  shifts  in  (i)  the  rate  at  which  CsCl  was  supplied 
to  the  flame,  (ii)  throughput  of  flame  gases  to  the  mass  spectrometer  or 
(iii)  mass  spectrometer  sensitivity  all  may  occur  during  a run,  large  devia- 
tions occasionally  occurred  in  the  set  of  values  for  kj  from  a single  ion  pro- 
file, in  which  case  the  runs  were  repeated.  Figure  1 is  an  Arrhenius  plot  of 
the  average  values  of  ki  found  for  each  flame.  Each  point  represents  the 
average  of  values  obtained  from  four  to  six  ion  profiles.  The  experimental 
uncertainty  in  each  point  is  % i 3.  The  least  squares  line  fitted  to  the 
data  gives 

k,  = 3 x i<r10  exp(-13100/T)  ml  molecule"1  sec"1  (6) 

The  question  of  the  validity  of  the  simplifying  assumption  used  in  writing 
Eq.  (4)  may  be  answered  in  part  by  the  following  arguments.  With  regard  to 
diffusion,  rough  calculations  show  axial  diffusion  to  have  a very  small  effect 
on  the  axial  distribution  of  species  in  the  flame.  A 'diffusion  velocity', 
v^  may  be  defined  by 


vd 


D d [Cl~] 

[Cl-]  dz 


(7) 


where  D is  the  ambipolar  diffusion  constant  of  Cl-.  The  maximum  value  of 
[Cl-]-1  (d [Cl-] /dz)  observed  is  % 2 cm-1.  D^-  in  these  flames  may  be  esti- 

mated to  be  ^ 10  cm2  sec-1  (Dh3q+  12  cm2  sec-1  in  flames  such  as  these10). 
Thus  the  ratio  v^/v^  ^ 0.02,  i.e.,  axial  diffusion  would  require  no  more  than 
2%  correction  to  the  data.  Similar  considerations  lead  to  similar  results 
with  regard  to  radial  diffusion. 
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The  possibility  that  the  reaction 

e + Cl  + M -*•  Cl-  + M 

may  be  competing  with  Reaction  (1)  should  be  considered.  The  ratio  of  the 
rate  of  Reaction  (8)  to  that  of  Reaction  (1)  would  be  given  by  the  expression 

k8  [Cl]  [M]  * ka  [HUM] 

kTlici]  kTKl0“l hT  (9) 

where  K10  is  the  equilibrium  constant  for  the  reaction 

H + HC1  t H2  + Cl  (10) 

which,  because  of  the  large  concentration  of  H and  H2 , must  be  nearly  balanced 
throughout  the  burned  gas  region  of  the  flames.  Values  for  K! 0 [H] [M]H2 ]-1 
(in  order  of  decreasing  magnitude) are  listed  in  Table  II.  From  this  data  it 
may  be  seen  that  the  contribution  to  the  rate  of  production  of  Cl-  from 
Reaction  (8)  would  be  greatest  in  the  760  Torr,  2475  K and  2250  K flames  and 
the  100  Torr,  1980  K flame.  These  flames,  then,  should  yield  anomalously  high 
values  for  kv  relative  to  the  other  flames  if  Reaction  (8)  was  important.  The 
two  flames  which  gave  values  of  k2  furthest  from  the  least  squares  line  in 
Fig.  1 are  the  760  Torr,  2250  K flame  (low)  and  the  760  Torr,  1818  K flame 
(high) . This  is  opposite  to  what  would  be  expected  if  Reaction  (8)  were  con- 
tributing to  Cl”  formation.  In  fact,  the  data  of  Fig.  1 show  no  tendency  at 
all  to  be  sensitive  to  the  parameter  in  the  right  hand  column  of  Table  II. 

Thus  the  data  appear  to  argue  against  Reaction  (8)  being  a significant  process. 

Recent  cross  section  measurements3’12  for  Reaction  (1)  show  that  the  peak 
cross  section  occurs  at  81  kJ  mole-1  (9.84  eV) . On  the  basis  of  these  deter- 
minations the  activation  energy  for  Reaction  (1)  measured  here  (109  kJ  mole-1) 
appears  too  large.  However,  it  may  be  seen  from  the  dashed  line  in  Fig.  1 
which  corresponds  to  an  expression  for  k2 , with  an  activation  energy  of 
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81  kJ  mole-1  that  this  temperature  dependence  falls  well  within  the  limits  of 
error  imposed  by  the  scatter  in  our  data.  If  81  kJ  mole-1  is  the  actual  activa- 
tion energy,  the  expression  for  kt  becomes 

. ki  = 5 * 10-11  exp  ml  molecule-1  sec-1  (11) 

The  comparison  of  the  rate  coefficient  measured  here  with  values  inferred 
from  electron  beam  and  swarm  experiments  requires  evaluation  of  the  following 
integral 


k 


23/2 

/TT  , 3 3 v 1 / 2 

(Hmk  T ) ' 


Q E exp 


(-  m) 


dE 


(12) 


Where  Q is  the  cross-section,  m the  reduced  mass  (approx,  the  mass  of  the 
electron)  and  E the  center-of-mass  kinetic  energy.  The  shape  of  the  cross- 
section  vs.  energy  curve  peaks  sharply  at  E0  = 0.84  eV.3»12  The  width  of  the 
cross-section  is  not  known,  but  one  would  expect  onset  to  be  near  0.82  eV,  the 
value  of  the  difference  between  dissociation  energy  of  HC1  and  the  electron 
affinity  of  Cl.13  The  assumption  that  the  cross  section  can  be  represented 
by  a Gaussian 


q - Q. 

\ 2a2 


where  a = 0.02  eV  leads  to  the  following  estimation  of  ki  at  2000  K 


ki 


'V 


4E0o  Qq 

ml/2  (k3X3)*/a 


exp(-  ff)  = 


4*107  Qc 


molecule-1  sec" 


(13) 

If  Qo  is  taken  to  be  4 x 10-1B  cm2  or  2 * 10-17  cm2,  the  low1  and  high3  values 
for  the  measured  peak  cross  section,  pre-exponentials  of  8 x 10-1°  and 
1.6  x 10-1°,  respectively,  are  obtained  in  Eq.  (13).  Thus  the 
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present  measurement  appears  to  favor  somewhat  the  former,  lower  value,  being 
at  2000  K about  a factor  of  3 below  that  inferred  from  Ref.  1. 

Regardless  of  which  of  the  above  expressions  for  ki  is  adopted,  a consider- 
able discrepancy  exists  between  these  direct  high  energy  (temperature)  measure- 
ments of  the  rate  of  Reaction  (1)  and  the  values  derived  via  the  equilibrium 
constant  and  measurements  of  the  rate  coefficient  for  the  reverse  associative 
detachment.  For  example,  the  measurements  of  Burdett  and  Hayhurst5  and  Howard 
et  al*  of  k_i  at  lower  temperatures  lead  to  values  of  the  pre-exponential  for 
ki  of  4 * 10-7  to  2 x 10-6  when  extrapolated  to  2500  K.  As  noted  in  Ref.  5, 
these  values  indicate,  in  turn,  peak  cross  sections  for  Reaction(l)  on  the 
order  of  10-1*  cm2  or  larger — some  103  to  106  times  greater  than  those  of  the 
present  work  or  of  Refs.  1-3.  We  can  offer  no  explanation  for  this  discrepancy 
at  this  time  other  than  marked  changes  in  the  detailed  reaction  mechanism  with 
temperature  or  a temperature  dependence,  for  the  pre-exponential  of  kt  such  as 
T-3.  These  questions  are  likely  to  be  resolved  only  by  direct  measurements 
of  ki  at  intermediate  temperatures. 
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TABLE  I.  FLAME  CHARACTERISTICS 
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Fig.  1 Average  values  for  the  dissociative  attachment  rate  coefficient 
obtained  in  the  seven  test  flames.  The  solid  line  is  a least- 
squares  fit  to  the  data  representing  ki  = 3 x 10-l°  exp(-13100/T) . 
The  dashed  line  has  a slope  of  81  kJ  mole-1,  the  activation  energy 
expected  from  the  work  of  Refs.  3 and  12. 
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This  paper  describes  briefly  a number  of 
plasma  diagnostic  techniques  which  have  been  widely 
used  for  the  measurement  and  characterization  of 
charged  species  in  combustion  systems.  The  limits 
of  applicability,  sensitivity  and  accuracy  of  the 
methods  are  discussed  and  the  working  equations  by 
which  the  experimentalist  reduces  data  are  pre- 
sented. Extensive  literature  references  are  also 
provided  in  which  to  find  detailed  accounts  of 
apparatus  construction  and  the  development  of  the 
mathematical  models  upon  which  the  working  equa- 
tions are  based.  Techniques  described  include 
electrostatic  probes,  microwave  cavity  and  cyclo- 
tron resonance,  electromagnetic-wave  attenuation 
and  mass  spectrometry. 

Introduction 

The  ionization  phenomena  associated  with 
nearly  all  practical  combustion  systems  have  been 
subjects  of  scientific  research  for  well  over  a 
century.  Subsequent  to  World  War  II,  appreciation 
of  the  importance  of  ion  and  electron  chemistry  in 
systems  such  as  rocket  exhausts,  re-entry  wakes, 
combustion-driven  HUD  devices  and  flame  Ionization 
detectors  lias  stimulated  the  development  and  appli- 
cation of  a wide  variety  of  diagnostic  tools  for 
the  detection  and  measurement  of  charged  species  in 
combustion  environments.  The  objectives  of  this 
paper  are  to  review  a number  of  the  more  commonly 
employed  diagnostic  techniques,  illustrate  their 
utility  by  presenting  samples  of  experimental  data 
and  provide  potential  users  with  literature  refer- 
ences to  more  detailed  descriptions  of  their  appli- 
cations. The  methods  discussed  in  greatest  detail 
are  those  which  have  proven  most  useful  to  the 
laboratory  experimentalist  as  tools  rather  than  as 
subjects  of  research  in  themselves.  It  is  not  the 
intention  of  this  paper  to  present  an  extensive 
description  of  the  theoretical  development  under- 
lying various  methods;  that  information  is  included 
only  by  reference  to  the  literature.  The  emphasis 
here  is  on  the  capabilities,  limitations  and  accu- 
racy of  a few  well-used  techniques,  and  the  working 
equations  through  which  the  raw  data  obtained  are 
converted  to  the  desired  appropriate  plasma  proper- 
ty. The  discussion  includes  descriptions  of  elec- 
trostatic probes,  some  microwave  and  radiofrequency 
methods  and  mass  spectrometry;  the  various  optical 
techniques  used  for  combustion  plasma  studies  (i.e. 
methods  employing  wavelengths  up  to  and  including 
the  far  Infrared)  are  not  discussed  here;  the 
Interested  reader  is  referred  to  the  spectroscopic 
literature  (e.g.  Rcfs.l  and  2). 


upon  the  magnitude  of  the  voltage  differential  and 
the  relative  sizes  of  the  probes  and  the  combustion 
system.  If  the  impressed  electric  field  is  large 
(-  100  V/cm) , and  the  probes  are  substantial  in 
size  (i.e.  approaching  the  size  of  the  flame,  or 
larger),  charged  species  will  be  drawn  out  of  flame 
essentially  as  they  are  formed.  Recombination  is 
thus  all  but  eliminated  and  the  current  measured 
between  the  electrodes  is  directly  proportional  to 
the  ion  formation  rate.  These  measurements  serve 
as  the  basis  for  the  flame  ionization  detectors 
used  primarily  in  gas  chromatography  for  the  detec- 
tion and  identification  of  hydrocarbons . A schema- 
tic of  a typical  detector*  of  this  kind  is  given  in 
Fig.  1.  It  depends  for  its  operation  on  the  low 
level  of  ionization  exhibited  by  Hj/air  flames*  and 
the  fact  that  the  combustion  of  non-ox  id i zed, 
carbonaceous  material  introduced  into  these  flames 
gives  rise  to  a nearly  constant  fraction*  of  ions 
per  carbon  atom  regardless  of  the  che.ilcal  struc- 
ture of  the  molecule.  These  systems  arc  routinely 
made  to  respond  linearly  to  the  carbon  content  of 
the  gas  stream  over  a 107-fold  range  of  sensitivity 
and  they  are  thus  ideally  suited  to  the  analysis  of 
gas  streams  emerging  from  chromatography  columns. 
The  probe  response,  however,  is  simply  a macro- 
scopic measure  of  flame  conductivity;  it  reveals 
nothing  about  the  details  of  the  gas  dynamics  or 
the  chemistry  of  the  flame. 
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Electrostatic  Probes 

The  oldest  and,  still,  the  most  commonly  per- 
formed flame  plasma  measurements  involve  a deter- 
mination of  the  current  flow  between  two  electrodes 
placed  in  or  near  the  flame  and  biased  at  different 
potentials.  Several  kinds  of  information  may  be 
extracted  from  data  obtained  in  this  way,  depending 


Figure  1.  Cross-Sectional  View  of  Flame  Ioniza- 
tion Detector  for  Gas  Chromatography, 
(from  Ref.  4) 

The  use  of  very  small  probes  at  variable 
collection  voltages  has  produced  the  greatest 
proportion  of  out  present  knowledge  of  the  detailed 
plasma  properties  of  high  temperature  combustion 
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systems.  The  current  vs.  voltage  ’probe  curves' 
obtained7  from  small  cylindrical  or  spherical 
probes  can,  in  principle,  be  interpreted  to  yield 
local  values  of  charged  species  concentrations  and 
electron  temperatures.  The  probes  themselves  are 
fairly  easy  to  build;  the  minimum  necessary  instru- 
mentation includes  only  a variable  d.c.  voltage 
source  and  an  electrometer  of  some  sort.  The  data 
is  deceptively  simple  to  collect  and,  if  anything, 
there  is  an  overabundance  of  theories  available 
with  which  to  interpret  the  results. 

Double  probes,  as  the  name  implies,  comprise 
two  electrodes  of  the  same  size  inserted  into  the 
plasma  for  the  purpose  of  withdrawing  charged 
species.  Since  the  critical  measurement  is  the 
flow  per  unit  probe  area  of  ions  or  electrons 
across  the  gaseous  sheath  of  shifting  potential 
surrounding  the  probe,  it  is  extremely  Important 
that  voltage  losses  from  the  probe  through  its 
insulation  be  eliminated  and  that  the  probe  area  be 
well  defined.  The  construction  of  such  probes  for 
use  in  hostile  combustion  environments  is  not  a 
trivial  task. 

Figure  2 illustrates  a typical7  probe  of 
cylindrical  geometry.  The  quartz  tubing  provides 
the  insulation  which  is  protected  from  dielectric 
degradation  due  to  heating  by  the  silver  cooling 
Jacket.  The  probe  wire  itself  must  be  of  a mate- 
rial refractory  to  flame  constituents  and  capable 
of  heating  to  high  temperatures.  It  should  also  be 
non-catal ytic  with  respect  to  atom  and  radical 
recombination  to  preclude  overheating  due  to  these 
effects.  Platinum*  and  (prefetably)  platinum-rhodium 
alloys  (e.g.  see  Kefs.  7,  9-11)  have  been  used  suc- 
cessfully by  a number  of  workers;  it  is  vise  to 
avoid  the  use  of  platinum,  however,  since  it  is 
notorious  for  exhibiting  catalytic  effects  in  H- 
containing  flames.  It  should  also  be  well  noted 
that  the  probes,  whether  'double*  or  'single',  must 
be  kept  as  small  as  practicable  so  that  the  probe 
operation  itself  does  not  alter  the  plasma  proper- 
ties it  seeks  to  measure;  i.e.  the  probe  should  not 
withdraw  enough  charged  species  to  change  the 
plasma  potential  significantly. 


Figure  2.  Electrically  and  Thermally  Insulated 
Probe. 

An  example  of  a cur rent -voltage  curve  from  a 
double  probe  is  given  in  Fig*  3.  It  can  be  seen 
from  this  data  that  the  'saturation'  currents  at 
high  potentials  dr^wn  by  the  two  probes  are  nearly 
equal  and  that  the  curve  is  nearly  symnetr ical . 

This  is  due  to  the  fact  that  the  charged  species  of 
lowest  mobility  (greatest  mss)  controls  trio  cur- 
rent flow  to  both  probes  and  only  its  conrentra t ion 
nay,  therefore,  be  measured.  In  high  temperature 


systems  then,  where  the  predominant  negatively 
charged  species  are  electrons,  these  curves  are 
controlled  by  positive  ion  collection  and  the 
method  is  Incapable  of  yielding  electron  densities. 
A perfectly  symmetrical  double  probe  curve,  un- 
fortunately, represents  an  ideal  situation  in  which 
the  probes  are  identical  in  size  and  geometry  and 
both  are  surrounded  by  gases  with  identical  plasma 
potentials.  In  practice  (e.g.  in  the  case  illus- 
trated by  Fig,  3),  this  condition  is  difficult  to 
achieve  for  a variety  of  reasons  discussed  in 
detail  in  Ref.  9.  Despite  these  difficulties, 
double  probe  methods  have  been  widely  and  success- 
fully employed  to  measure  positive  ion  concentra- 
tions and,  from  the  nature  of  the  curve  in  the 
rapidly  changing  portion  at  email  probe  voltage 
differences,  electron  temperatures.  The  means  by 
which  the  data  is  interpreted  in  both  these  cases, 
however,  is  identical,  in  essence,  to  that  employed 
for  single  probes  and  they  can  be  conveniently 
discussed  jointly. 


Figure  3.  Double  Probe  Current-Voltage  Character- 
istic (from  Ref . 8) . 

If  the  size  of  the  positively  charged  electrode 
can  be  increased  sufficiently,  saturation  with 
respect  to  the  electron  current  can  be  achieved  and 
electron  densities  can  be  derived  as  well.  This 
condition  is  attained  by  making  the  second  electrode 
the  burner  and/or  large  screens  inserted  at  various 
positions  in  the  plasma,  taking  care  that  conditions 
about  the  small  probe  arc  not  affected.  The  ratio 
of  the  probe  arras  required  to  achieve  adequate 
positive  ion  transport  when  the  small  probe  is^ 
sampling  electrons  has  been  shown  3 to  be  - 10  . 

This  is  equivalent  to  saying  that  the  'single* 
probe  must  be  very  Mnall.  The  difficulty  in 
achieving  this  condition  Is  illustrated  in  Fig.  4 
which  is  a sample  single  probe  curve  of  a kind 
unfortunately,  frequently  encountered;  electron 
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Figure  4.  A Typical  Single-Probe  Current-Voltage 
Characteristic  (from  Ref.  9). 

saturation  at  positive  probe  potentials  has  not 
been  attained  due  to  the  inadequate  size  of  the 
negative  electrode.  The  large  electrode  in  this 
case  was  a 3.8  cm  diameter  Mekcr-type  burner, 
giving  an  electrode  area  ratio  of  '■'*  150.  An 
alternative  technique  to  increase  the  positive 
ion  flux,  involving  the  addition  of  a third 
independently  biased  electrode,  has  been  suggested 
by  Yanagc*3’  lh  but  there  are  reservations  2 as 
to  its  effect  on  the  bulk  plasma  properties.  In 
short,  although  electron  concentration  measurements 
using  single  probes  are  straightforward  in  princi- 
ple, they  are  very  difficult  in  practice  and 
great  care  must  be  exercised  in  order  to  obtain 
accurate  results. 

Many  theories  have  so  far  been  advanced  for 
the  interpretation  of  electrostatic  probe  curves 
taken  under  various  conditions.  The  theories 
take  into  account,  anong  many  things,  probe 
geometry,  the  magnitude  of  the  applied  potential, 
the  effects  of  plasma  velocity  and  the  changes  in 
plasma  properties  induced  by  variations  in  combus- 
tion system  pressure.  To  the  laboratory  experi- 
mentalist this  frequently  represents  an  overabun- 
dance of  models  from  which  to  chouse;  in  fact, 
the  differences  in  charged  species  concentrations 
computed  using  the  various  approaches  seldom  fall 
outside  the  accuracy  limitations  imposed  by  the 
data  itself.  One  of  the  simplest  and  most  widely 
used  theories  for  the  interpretation  of  probe 
data  In  combustion  systems  is  that  of  Bohn, 

Burhop  and  Massey1'  as  modified  by  Calcote. 1 * * 1 7 
This  theory, for  the  usual  case  of  cylindrical 


3 


probes  (of  length  L and  diameter  d) , relates  the 
value  of  the  positive  ion  current  (j+)  obtained 

by  extrapolating  the  negative  voltage  portion  of 
the  probe  curve  to  the  plasm3  potential  (V  , see 

P 

Fig.  4),  to  the  positive  ion  concentration,  N , 
through  the  expression  + 


K( 


2 n m+  > 
k T+  ea> 


+ 


0.19  Ld 


where  ra+  is  the  positive  ion  molecular  weight,  T+, 
Its  temperature  and  X+,  its  mean-free  path;  e is 
the  charge  of  an  electron,  X - L + 2A+  and  B «= 
(Xa-(d+2X+)*}  l/  2 . For  probes  of  other  than 

cylindrical  geometry,  only  the  quantity  in  brackets 
changes  form;  a corresponding  expression  for  el- 
lipsoidal (or,  in  one  limit,  spherical)  probes  is 
given  in  Ref.  7.  Directly  analogous  expressions 
may  be  written  for  electron  concentration  (N  ) 
determinations  in  which  all  the  subscripted  quan- 
tities for  positive  ions  are  replaced  by  their 
corresponding  electron  values. 


The  uncertainties  inherent  in  the  use  of  s 

Eq.  (1)  lie  in  tke  local  identification  of  m , 
the  estimation  of  X+  and  in  extrapolation  oi  the 

ion  current  curve  to  V . The  latter  is,  of  J 

p , 

course,  dependent  upon  proper  identification  of  j 

the  'break*  point  in  the  electron  current  vs. 
voltage  curve  which  defines  the  plasma  potential. 

Obviously,  if  saturation  is  not  obtained,  as  in  ; 

Fig.  4,  the  proper  value  of  cannot  be  deter- 
mined. Fortunately,  since  the  slope  of  the  5 

positive  ion  saturation  curve  is  small,  this 

source  of  error  is  usually  negligible  in  measure-  \ 

ments  of  positive  ion  concentraions.  It  should  i 

he  noted,  however,  that  this  is  true  only  of  the 
positive  ion  current;  the  extent  to  which  the 

probe  falls  short  of  achieving  electron  saturation  j 

is  directly  proportional  to  the  resultant  error  • < 

in  Ne.  I 

Charged  species  transport  due  to  convective  . 

flow  is  not  accounted  for  in  Eq.  (1)  although 
such  effects  are  commonly  encountered  in  combustion 
systems.  Neglect  of  this  effect  can  be  of  consid- 
erably greater  import  than  the  errors  in  positive 
ion  determinations  so  far  discussed.  Jensen  and 
Kurzius18  have  presented  a modified  version  of  1 

Eq.  (1)  which  includes  a consideration  of  mass 
flow  and  is  generally  applicable  to  any  laminar, 
high  velocity,  subsonic  system.  For  the  frequently 
encountered  instance  of  a cylindrical  probe  with 
its  axis  normal  to  the  flow,  n+  is  given  by 


2 n m. 


M<TT) 


(Sh)  D+ 


(2) 


where  D+  is  the  ambipolar  dlffusivity  and  (Sh)  is 

the  perimeter-mean  Sherwood  number  given,  in 
terms  of  the  Reynolds  (Re)  and  Schmidt  (Sc) 
numbers,  by  the  expression 

(Sh)  “ 0.55  (Re),/7  - 2 ln{2(I.  + 2A+)/d 

+ 2A+))  O) 
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Teats  of  this  approach  have  so  far  been  limited 
to  positive  ion  determinations;  the  validity  of 
the  analogous  expression  for  n has  not  been 
Investigated.  e 

In  addition  to  yielding  electron  and  positive 
ion  number  densities,  electrostatic  probes  are 
also  used  for  measurements  of  electron  tempera- 
tures. However,  the  proper  means  of  performing 
these  measurements  has  been  the  subject  of  consider- 
able controversy  for  more  than  a decade.  In 
principle,  the  slope  of  a plot  of  ln[electron 
current]  vs.  applied  voltage  is  equal  to  e/kT  , 
the  electron  current  being  the  observed  external 
probe  circuit  current  minus  the  extrapolated 
positive  ion  saturation  current  at  the  applied 
voltage.  Workers7*19*20  using  single  probe 
curves  in  this  fashion,  however,  frequently  found 
extremely  high  electron  temperatures  in  a narrow 
region  within  hydrocarbon  flames  while  those  employ- 


ing double  probes  * 


in  similar  flames,  usually 


obtained  values  close  to  the  gas  temperature 
throughout  the  entire  flame.  It  is  now  generally 
accepted  that  the  latter  is  true;  the  departures 

of  from  the  gas  temperature  in  normal  f Limes 

are,  at  most,  a few  hundred  degrees* * » 2 3 and  they 
only  occur  very  close  to  the  flame  front.  The 
source  of  these  discrepancies  in  single  probe 
measurements  are  attributed  to  the  use  of  data 
taken  too  far  from  the  positive  ion  saturation 
region;  when  single  and  double  probe  data  from  the 
portions  of  the  respective  curves  close  to  positive 
ion  saturation  are  directly  compared,  as  in  Ref. 

9,  the  T^’s  are  nearly  identical.  Similarly, 

despite  the  diverse  assumptions  underlying  the 
various  available  probe  theories  and  the  supposed 
differences  in  their  domains  of  applicability,  the 
same  conclusion  is  reached24*9  using  the  simple 
theory  of  Calcote7  or  those  of  Su  and  Lam,27 
Lam,2*  Johnson  and  Walter,29  or  Cohen.20  For  the 
purist  who  wishes  to  analyze  electrostatic  probe 
behavior  in  detail,  extensive  discussions  of  the 
above,  and  many  other  theories  are  conveniently 
brought  together  in  Ref.  31. 
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Figure  5.  Burner  and  Microwave  Cavity. 

A representative  circuit*2  for  utilizing 
cavity  resonance  effects  is  given  in  Fig.  6. 
Microwave  power  from  a tunable  klystron  (V  ) with 
a variable  output  impedance,  Rg,  Is  fed  in^:o  the 

cavity  through  one  of  the  antennae  indicated  in 
Fig.  5.  L,  C,  and  represent,  respectively,  the 

inductance,  capacitance  and  resistance  of  the 
cavity  itself . The  impedance  is  provided  by  a 

crystal  diode  rectifier  across  which  a d.c.  voltage 
(V^)  or  current  flow  (I)  is  measured.  After 

setting  up  the  circuit  and  determining  the  TM(OXO) 
resonant  frequency  in  the  presence  of  a flame 
containing  no  electrons,  an  ionizable  material  may 
be  added  to  the  flame  to  produce  electrons  and  the 
characteristics  of  the  circuit  redetermined.  An 
analysis  of  the  electrical  properties  of  the 
cavity  shows  that  the  Q of  the  cavity  and  its 
resonant  frequency  before  the  addition  of  electrons 
(subscript  'o')  is  related  to  that  with  electrons 
(no  subscript) , through  the  expression 


Microwave  Cavity  Resonance 

If  a microwave  cavity  is  built  so  as  to  in- 
clude a plasma  or  a flame  containing  ionized 
species,  its  Q factor  (defined  as  the  ratio  of 
energy  stored/energy  lost)  at  resonance  becomes  a 
function  of  the  dielectric  constant,  and  thus  the 
electron  content,  of  the  plasma.  Specifically, 
the  technique  which  has  proven  most  useful  is  one 
in  which  the  cavity  is  of  a geometry  suitable  for 
resonance  in  the  TM(010)  mode  at  a specific  micro- 
wave frequency,  u,  between  l and  4 GHz.  A sche- 
matic of  such  a cavity,  together  with  the  burner 
and  flame  system  to  which  it  is  coupled,  is  shown 
in  Fig.  5.  In  practice, 3 2-35  these  cavities  are 
annuli  about  10  cm  In  diameter  and  a few  mm  high. 
The  hole  in  the  center  can  be  as  large  as  about  3 
cm  without  seriously  affecting  it3  performance. 

They  must,  of  course,  be  cooled  to  prevent  thermal 
degradation  by  the  flame  and  they  should  be 
therroostated  to  minimize  drift  due  to  small  changes 
in  geometry.  The  spatial  resolution  of  electron 
concentration  measurements  made  along  the  flame 
axis  is  typically  about  the  same  as  the  cavity 
height  except  in  the  region  of  the  flame  close  to 
the  burner  where  the  Q factor  of  the  cavity  is 
Influenced  by  the  presence  of  the  burner  Itself. 


fr-1 


where  A is  a constant  dependent  on  cavity  geometry 
and  o is  the  conductivity  of  the  plasma.  For 
caviti.es  of  the  sort  described  above,  it  has  been 
shown  6 shown  that  the  frequency  shift,  in  reality, 
is  small  in  comparison  to  the  change  in  o and  the 
term  in  brackets  in  Eq.  (4)  can  be  neglected. 

Now,  if  Rgfcj/CRg  + R^  >>  Rc»  it  to*5  beco  further 

shown  that  I « Q2  and,  since  a « n , one  may 
derive  the  following  expression  for  the  electron 
concentration  in  terms  of  the  d.c.  currents  before 
and  after  their  introduction  (I0  and  1, 
respectively) : 


In  this  equation,  C is  a number  of  combined  con- 
stants and  must  be  determined  by  direct 
calibration. 
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Figure  6.  Equivalent  Circuit  of  Cavity  System 
(from  Ref.  32) . 

To  calibrate  the  cavity,  a flame  of  known 
electron  content  can  be  used.  A large  number  of 
these  are  now  available  as  a result  of  an  extensive 
series  of  investigations*  using  cavities  and  a 
number  of  other  diagnostic  methods  as  comparative 
cross-checks.  Once  a value  of  C has  been  deter- 
mined from  the  known  n and  the  observed  d.c. 
e 

current  ratio,  it  is  valid  over  the  entire  avail- 
able range  of  cavity  response  but,  strictly 
speaking,  only  for  that  flame  used  in  the  calibra- 
tion. Since  C depends  on  the  electron-neutral 
collision  frequency  ( vc) , which  is  a function  of 

gas  composition  and  temperature,  it  will  vary  from 
flame  to  flame.  This  is,  fortunately,  not  found, 
to  be  a significant  source  of  error  unless  the  H20 
content  of  the  flame  changes  appreciably;  a 
calibration  performed  in  H*/02  flames  would  not, 
however,  be  expected  to  be  adequate  in  flames  of 
C0/0a,  for  example. 

The  accuracy  of  this  technique  is  as  good  as 
the  calibration  allows  and  the  sensitivity  is  a 
direct  function  of  the  cavity  constant.  In  its 
implementation,  the  method  has  proven  useful, 
practical ly,  in  the  range  10°  < n < 10* 2 cra~3. 

At  the  * lower  limit  of  this  range,  the  system  must 
be  very  noise-free  and  stable;  even  with  well- 
designed  cavities,  at  n --108  cri"3,  one  must  be 
able  to  detect  a change  in  I of  a few  parts  per 
103  or  less.  The  most  severe  limitations  on  the 
usefulness  of  the  method  are  imposed  by  its  geomet- 
rical requirements.  The  plasma  to  be  studied  must 
be  a part  of  the  cavity  and  thus  the  overall 
dimensions  of  the  apparatus  must  be  large  relative 
to  those  of  the  flame.  Despite  its  applicability 
to  rather  small  combustion  systems  of  somewhat 
specialized  conf iguration,  this  technique  lias 
proven  of  great  value  in  basic  laboratory  combus- 
tion studies. 

Cyclotron  Resonance 


When  a flame  or  other  plasma  is  intersected 
by  a magnetic  field  and  a microwave  beam  at  right 
angles  to  one  another,  absorption  cf  electrical 
energy  may  occur  and  the  transmitted  electrical 
energy  is  attenuated.  In  low  pressure  plasmas, 
in  which  the  microwave  frequency,  u,  Is  made  large 
with  respect  to  the  plasma  frequency,  the 

maximum  attenuation  takes  place  at  the  electron 
cyclotron  resonant  frequency.  The  nagnitudr  of 
this  attenuation,  and  the  bandwidth  over  which  it 
occurs,  arc  directly  related  to,  and  thus  can  be 


used  to  measure,  both  ng  and  the  electron  collision 
frequency,  vc. 

Since  the  impression  of  a significant  magnetic 
field  on  the  plasma  requires  rather  bulky  apparatus, 
this  method  is  also  applicable  primarily  to  small 
laboratory  systems.  The  most  convenient  measure- 
ments have  proven*0 "**  to  be  those  in  which  the 
beam  frequency,  w,  is  fixed  at  a value  >>  w , the 
magnetic  field  strength,  H,  is  swept  througt!  the 
resonant  condition  and  the  microwave  power  trans- 
mission is  measured  as  a function  of  H.  A plot  of 
0,  the  microwave  attenuation,  vs.  H then  exhibits 
a peak  at  some  field  value,  H0,  which,  together 
with  the  measured  width  of  the  peak  at  half  maximum, 
gives  the  electron  collision  frequency  via  Eq.  (6): 


v 

e 


AH  to 
2 Ho 


(6) 


For  this  equation  to  be  valid,  the  width  of  the 
resonance  line,  AH,  must  be  the  result  of  colli- 
sional  broadening  and  this  requirement  is  very 
stringent.  Gray*5  has  computed  that,  in  fact, 
(u>p/id)2  must  be  'considerably  below  0.1'.  This 

restriction,  in  practice,  limits  the  applicability 
of  the  method  to  subatmospher ic  pressure  plasmas 
and  flames  with  1 'v  1011  cnT3. 

Within  the  same  limitations  as  above,  the  B 
vs.  H curve  yields  the  electron  concentration 
through  the  equation*6 

n = 3.8  x 107  0 AH/L  (7) 

e 

where  L is  the  microwave  path  length  through  the 
(uniform)  plasma.  This  measurement,  and  the 

determination  above,  when  made  via  this  technique, 
have  the  advantages  that  no  perturbation  of  the 
plasma  by  a probe  is  necessary  and  that  the  method 
requires  no  calibration.  Its  accuracy,  in  princi- 
ple, should  be  quite  good  since  no  absolute  mea- 
surements are  required  other  than  those  of  u>  and 
H,  both  of  which  can  ho  quite  precisely  metered. 
Resides  the  rather  limited  scope  of  applicability, 
however,  the  method  also  suffers  the  disadvantage 
of  poor  spatial  resolution  (y  0.5  cm)  so  that 
detailed  flame  structure  determinations  are 
difficult . 

Direct  Attenunt ion  and  Other  Microwave  Methods 

Since  both  of  the  resonance  methods  described 
above  involve  surrounding  a combustion  system  with 
a considerable  amount  of  apparatus,  the  sizes  of 
the  plasmas  that  are  accessible  to  study  is  limited, 
in  reality,  to  those  of  laboratory  test  flames. 

The  plasnu)  properties  of  larger  systems  such  as 
rocket  exhausts  must  be  examined  through  the  use 
of  rugged  electrostatic  probes*7* *8  or  through  the 
interaction  of  the  plasma  with  focused  beams  of 
electromagnetic  energy,  specifically,  riiorovivcs. 
When  a microwave  beam  encounters  a plasma,  it  can 
be  measurably  attenuated,  refracted  or  reflected, 
depending  on  the  local  values  of  n^,  and  w. 


* The  series  includes  work  by  a number  of  workers  and  is,  in  facL  so  extensive  that  it  is  impractical  to 
list  them  here.  The  interested  reader  is  referred  to  ^eviews  of  the  subject,  e.g.  Refs.  37-39. 
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Attenuation  measurements  must  be  made  under 

the  conditions  u >>  w ; if.  in  addition,  u £ v the 
P ® 

attenuation,  0,  in  dB/cm  path  length  is  given  by 


Values  of  n and  v can,  therefore,  both  be  ob- 
e e ’ 

tained  from  a determination  of  attenuation  at  two 
or  more  different  frequencies.  A schematic  of  a 
typical  apparatus  for  measuring  attenuation  in  a 
rocket  exhaust  flame  is  shown** * in  Fig.  7. 
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Figure  7.  Apparatus  for  Attentuation  Measure- 
ments in  Rocket  Exhausts  (from  Ref.  49). 


Attenuation  measurements  have  been  used  to 
study  small  laboratory  test  flames30"54  as  well  as 
rocket  exhausts.  However,  when  applied  to  small 
systems,  wavelength  considerations  severely  limit 
the  spatial  resolution  possible.  In  addition,  an 
examination  of  the  frequency  requirements  outlined 
above  reveals  that,  regardless  of  the  system  under 
study,  the  attenuation  method  is  applicable  only 
when  ng  X 10*°  cm“*»  For  these  reasons,  cavity 
resonance  has  essentially  replaced  attenuation  in 
the  study  of  laboratory  scale  flames. 

Microwave  reflection  (scattering)  represents 
an  alternative  method  for  characterizing  the  elec- 
trical properties  of  large  scale  combustion 
systems. 33' 36  The  apparatus’7*38  for  these  mea- 
surements is  similar  in  design  to  those  employed 
Jn  attenuation  determinations  and  the  applicable 
pan&e  of  wavelengths  is  identical.  This  technique 
is  a difficult  one  to  apply,  however,  being  very 
sensitive  to  the  conditions  at  the  plasma  bound- 
aries.3**50 In  general,  it  is  sufficiently  accu- 
rate only  for  order  of  magnitude  estimates  of  J7 
at  concentrations  of  101J  cm*'3  and  higher.  Measure- 
ments of  electron  densities  via  determination  of 
the  refractive  index  are  comparably  insensitive**5 
and  somewhat  more  complex  in  the  execution.  They 
have  seen  very  limited  use. 

Mass  Spectrometry 

A large  proportion  of  the  most  recent  detailed 
mechanistic  studies  of  both  ionic  and  neutral 
species  processes  in  combustion  systems  have  been 
performed  using  mass  spectrometers  equipped  with 
specially  designed  sampling  systems.  These  instru- 
ments, whether  applied  to  studies  of  neutrals  or 
ions,  share  the  advantages  inherent  to  mass  spec- 
trometry and  the  problems  associated  with  obtaining 


representative,  chemically  frozen,  gas  samples.  A 
discussion  of  the  molecular  beam  sampling  systems 
used  to  investigate  neutral  species  is  beyond  the 
scope  of  this  paper  and  the  interested  reader  is 
instead  referred  to  some  examples  in  the  litera- 
ture.53-54 In  principle,  however,  the  experimental 
objectives  are  the  same:  one  must  withdraw  a gas 

sample  from  the  system  and  then,  somehow,  terminate 
chemical  reactions  and  separate  the  species  of 
interest  from  the  remainder  of  the  sample.  The 
first  of  these  goals  is  accomplished  by  expanding 
the  g as  through  a nozzle  or  orifice  into  a vacuum. 
The  reduction  in  pressure  accompanying  the  expan- 
sion affectively  freezes  all  but  the  most  rapid 
chemistry  by  decreasing  the  molecular  collision 
frequency;  the  corresponding  decrease  in  tempera- 
ture also  retards  all  chemical  reactions  except 
those  with  little  or  no  activation  energy.  In 
molecular  beam  work,  the  central  portion  of  the 
c-xpar.ding  gas  jet  (comprising  molecules  which  have 
not  undergone  collisions  with  the  walls)  is  then 
separated  from  the  remainder  of  the  gas  by  a 
'skimmer*  which  allows  the  molecules  near  the  axis 
to  flow  into  a second  chamber  containing  the  mass 
spectrometer  and  its  ion  source.  In  the  sampling 
of  ions,  the  function  of  the  skimmer  is  performed 
by  electrostatic  lenses  shown  scheratical ly55  in 
Fig.  8.  The  whole  process  of  ion  sampling  is  thus 
simplified  by  two  factors:  the  absence  of  gas 
dynamic  problems  associated  with  the  skimmer  and 
the  elimination  of  the  ion  source. 
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Figure  8.  Flame  Ion  Mass  Spectrometer 
(from  Ref . 65X 

The  sampling  orifice  in  typical  instruments 
is  usually  placed  at  the  t Ip  of  a cone.  If  con- 
densable materials  are  to  be  present  in  the  flame, 
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the  cone  is  constructed  of  a material  (e.g. 
quartz*1  or  stainless  steel,*1) which  will  not 
conduct  heat  away  rapidly  thus  allowing  the  tip  to 
become  hot  and  prevent  condensation.  The  cone 
angle  represents  a compromise  between  being  as 
large  as  possible  on  the  inside  to  facilitate 
pumping  and  ion  acceleration  and  as  small  as 
possible  on  the  outside  to  minimize  disturbance  to 
the  flow  of  flame  gases.  In  practice,  overall 
cone  angles  of  about  90*  are  typical.  The  extent 
to  which  the  flame  chemistry  is  frozen  is  the 
major  concern  in  sampling.  An  analysis  of  this 
problem  by  Hayhurst  and  Telford**  using  a simple 
gas  dynamic  model  of  the  sampling  process, 
indicates  the  critical  available  time  during  which 
unwanted  chemical  reactions  (i.e.  those  which 
occur  under  conditions  not  representative  of  the 
combustion  system,  e.g.  at  lower  temperature)  can 
occur  is  of  the  order  10“  sec.  This  number 
provides  a basis  for  estimating  whether  or  not 
significant  errors  of  this  nature  are  possible. 

For  example,  in  the  instance  of  an  ion-molecule 
reaction  between  an  ion  at  a concentration  of  10* 0 
cm-3  and  a neutral  free  radical  at  1015  cm“*  (both 
of  which  are  typical  of  flames  at  1 atm),  signifi- 
cant alteration  of  the  ion  concentration  during 
sampling  will  occur  if  the  rate  constant  for  the 
reaction  is  larger  than  about  10~8  era1  molec”1 
sec”1 . This  is  a large  value  not  commonly  en- 
countered even  among  ion-molecule  reactions.67*6* 

Once  ions  have  been  sampled  and  focused  into 
a mass  spectrometer,  the  nature  of  the  apparatus 
is  unimportant  except  that  its  throughput  efficiency 
will  determine  the  signal  strength  at  the  ion 
detector  and  the  system  being  studied  should  be 
shielded  from  any  strong  fields  associated  with 
its  operation.  The  pioneering  research  on  ioniza- 
tion in  flames  was  performed  with  magnetic*9  »70 
and  r.f.7  instruments  but  now,  quadrupole  mass 
filters  are  the  most  commonly  employed66*  61*  7 **  71 
because  of  their  high  transmission  efficiency  and 
readily  controllable  resolution  (one  is  always 
obtained  at  the  expense  of  the  other).  Both 
positive  and  negative  ions  are  routinely  studied 
in  these  mass  spectrometers  and,  by  matching  the 
initial  orifice  size  to  the  combustion  system 
pressure,  flames  at  a few  Torr  to  1 atm  pressure 
can  be  sampled. 

The  instrument  parameter  which  imposes  the 
greatest  limitations  on  overall  mass  spectrometer 
response  is  the  ion  transmission  through  the 
region  immediately  downstream  of  the  sampling 
orifice.  Most  of  the  ion  losses  occur  here  through 
collisions  with  the  walls  of  the  system.  Because 
of  the  non-reproducibility  of  these,  effects  and 
the  sensitivity  of  electrostatic  focusing  fields 
to  pump  oil  an  1 dirt  deposits  on  lens  components, 
absolute  calibration  of  the  apparatus  is  nearly 
impossible.  It  is  feasible  to  use  comparison  of 
an  ion  current  to  that  produced  by  a known  concen- 
tration of,  say,  alkali  metal  ions,  but  even  this 
procedure  must  be  done  carefully  to  account  for 
such  things  as  changes  in  mass  filter  throughput 
and  multiplier  response  with  mass.  Despite  these 
shortcomings  and  the  unavoidable  complexity  of  the 
apparatus,  mass  spectrometry  is  still  the  most 
powerful  and  versatile  tool  available  for  the 
study  of  corabusLion  plasma  chemistry. 
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